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Famous Words: 


Mathematics, rightly viewed, posses not only truth but supreme beauty — a beauty cold 
and austere, like that of sculpture. 


By Bertrand Russell, a British philosopher and mathematician 
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On a Boundary Value Problem with 


Fuzzy Forcing Function and Fuzzy Boundary Values 


Hiilya GULTEKIN CITIL 
(Department of Mathematics, Faculty of Arts and Sciences, Giresun University, Giresun, Turkey) 


E-mail: hulyagultekin55@hotmail.com, hulya.citil@giresun.edu.tr 


Abstract: In this study, a problem with fuzzy forcing function and fuzzy boundary values is 
investigated. The problem is solved by two different solution methods. Theorems are proved 
about solutions. Comparison results are given. Example is solved on studied problem. 
Graphics of the solutions are drawn. Conclusions are given. It is stated which method is 


more useful. 


Key Words: Fuzzy boundary value problems, second-order fuzzy differential equation, 


generalized differentiability. 
AMS(2010): 03E72, 34A07. 


§1. Introduction 


Fuzzy logic is studied by many researchers [10, 18]. In recent years, the topic of fuzzy differential 
equations has been rapidly growing [1, 11, 13, 20, 23]. Because, solving the fuzzy differential 
equations is a very important topic. Fuzzy differential equations can be studied by different 
approaches. These are Hukuhara differentiability [5, 15], generalized differentiability [2, 3] and 
to generate the fuzzy solution from the crips solution. There are at most four solutions for fuzzy 
boundary value problems using the generalized differentiability [17]. Liu [19] showed that these 
four solutions reduce to two different solutions when the function is monotone. To generate 
the fuzzy solution from the crips solution can be three ways. These are extension principle [5, 
6], the concept of differential inclusion [14] and fuzzy problem is to consider as a set of crips 
problem [8]. 

The aim of this study is to investigate the solutions of the fuzzy boundary value problem 


with fuzzy forcing function and fuzzy boundary values by two different solution methods. 


§2. Preliminaries 


Definition 2.1((22]) A fuzzy number is a mapping u: R — [0,1] with the following properties: 


(1) u is normal; 
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(2) u is convex fuzzy set; 
(3) wu is upper semi-continuous on R; 


(4) cl{x € R | u(x) > 0} is compact, where cl denotes the closure of a subset. 


Let Rp denote the space of fuzzy numbers. 
Definition 2.2((17]) Let ue Re. The a-level set of u is 
[ul ={2 ER| u(x) >a}, O<aK<l. 
[u]° = [u,,Ua] denotes the a-level set of u. 


Remark 2.1([7, 17]) The sufficient and necessary conditions for [u,,@q] to define the para- 


=a?) 


metric form of a fuzzy number as follows: 

(1) u, is bounded monotonic increasing (nondecreasing) left-continuous function on (0, 1] 
and right-continuous fora =0 , 

(2) %. is bounded monotonic decreasing (nonincreasing) left-continuous function on (0, 1] 
and right-continuous for a = 0, 

(3) u, <Ua,0<a<l. 


Definition 2.3([12, 17, 21]) Let u,v € Rp. If there exists w € Rr such that u=v+w, then 


w is called the Hukuhara difference of fuzzy numbers u and v,and it is denoted byw=uCv. 


Definition 2.4([4, 12, 17]) Let f : [a,b] > Rr and to © [a,b]. We say that f is Hukuhara 
differentiable at to, if there exists an element f (to) € Rr such that for allh > 0 sufficiently 
small, Sf (to +h) © f (to), f (to) O f (to — hk) and the limits hold 


lim 
h-0 h-0 


Flo SWS Fo) = tim LUO SF lfo — H) = f (to). 


Definition 2.5((17]) Let f : [a,b] — Rr and to € [a,b]. We say that f is (1)-differentiable at 
to, if there exists an element f (to) € Rr such that for all h > 0 sufficiently small (near to 0), 
exist f (to +h) © f (to), f (to) © f (to — h) and the limits 


am fo +4) OF lto) _  f (to) © F (to - h) 
h—0 h h—0 h 


=f (to), 


and f is (2)-differentiable if for allh > 0 sufficiently small (near to 0), exist f (to) Of (to +h), 
f (to —h) © f (to) and the limits 


f (to —h) © f (to) ' 


fig SOT Ce i — = f (to). 


h-0 —h h-0 


Theorem 2.6([16]) Let f : [a,b] + Rr be fuzzy function, where [f (t)|° = ie (t), f. (t)| for 
each a € [0,1]. 
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(i) If f is (1)-differentiable then f , and F,, are differentiable functions and f’ (| 
[fF], 

(ii) If f is (2)-differentiable then f and f.,, are differentiable functions and f On = 
fe Ost O|- 
Theorem 2.2({16]) Let f’ : [a,b] > Rr be fuzzy function, where [f (t)]* = li (t), fa (0) 
for each a € [0,1], f is (1)-differentiable or (2)-differentiable. 

(i) If f and f' are (1)-differentiable then f and F, are differentiable functions and 
[fo] =([£ 0.7.0), 

(«i) If f is (1)-differentiable and f' is (2)-differentiable then f' and F., are differentiable 
‘ol =Ro.£e], 


(iii) If f is (2)-differentiable and f' is (1)-differentiable then f and F, are differentiable 


functions and i 
functions and ba ()) = 7. (t) if (| , 

(iv) If f and f are (2)-differentiable then f. and q. are differentiable functions and 
7’) =([£0.7.0)- 


§3. Main Results 


Consider the two-point boundary value problem 


~N 


y (t)=A¥H+FE), yO)=8, WO =7, (3.1) 


are symmetric triangular fuzzy numbers and A > 0. 


3.1. Solution Method 1.((9]) Let divide the problem (3.1) into three different problems 


following: 


(i) The first problem is 


i O=wO+A, vO =225, wo = 2" (3.2) 


(ii) The second problem is 


y=), u0)= (S205), w= (G2054). ea) 
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(iii) The third problem is 


y =dy+(-1,0,1), y(0) =0, y(2 =0. (3.4) 


The solution of the differential equation in (3.2) is 


1 2 
y(t) = eet art = xt ~ 2° 


Using the boundary conditions, the coefficients c; and cz are found as 


d+d = é 
(4 + 30? 4 2) e vr (et? 4 2.) 
C= ’ 
eV — e- VE 
ye (cte 2 d+d , 192 2 
ea ($ #) (4 + x +2) 
C2 = ° 
eVX — e-VE 


Then, the solution of (3.2) is 


fie sinh (VA (€— 2) (at 2 


sinh (VX) ¥ ae 
jena (vx) dtd | 1 jp = 1 2 9 
CH) (st tos 8) ted 


Since 21 = eY*!, a» = e~Y™ are the linear independent solutions of the differential equation 
in (3.3), 


wi(t) = x (£) x1 (t) — x1 (£) x2 (t) 
1 (0) x2 (2) — x (£) x2 (0) 
sinh (vx (@— t)) 
7 sinh (v2) : 
wo (t) = 1 (0) i) (t) — 2x2 (0) rial (t) 
1 (0) x2 (2) — x (2) x2 (0) 
sinh (vxe) 
sinh (v22) 


Then, the solution of the problem (3.3) is 


von ay CFE) Seay CA) 
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Since the solution of the equation y = Ay-1lis 


4 sinh (VX(¢— 1) sinh (Vt) 
ee ‘{ Sinn (V2) smb (V0) 


and the solution of the equation y" =Ay+1is 


_ [sinh (Vr(e-0) _ sink (vx) ; 
i sinh (VX) ’ aa (vr) 


yi (t) = 


the solution of fuzzy boundary value problem (3.4) is 


y (t) = {min {y-1 (4) ,0, yr (¢)} ,0, max {ys (t) ,0, yn (é)}}- (3.7) 


Then, from (3.6) and (3.7), the fuzzy lower solution is 


sith (VAE-9) fez 1 | sinh (VM) fg 1\ 1 
uN sinh (v22) ( 2 ‘) sinh (v2.2) ( 2 : r’ on 


and the fuzzy upper solution is 


- sinh (VA (€— 2) ae _ sinh (VX) a ae ome 
7 (t) = sinh (V30) ( reas )4 ( es ) < (3.9) 


That is, the fuzzy solution is 


~N 


y (t) = (y(t),0,9(@). (3.10) 
Finally, from (3.5) and (3.10), the solution of the problem (3.1) is 


~N 


Y (t) = (y(¢), y(t). 9) (3.11) 


y(t) = aa) (c+ : as (a+ a4 + (0 :)) 


sinh (vr2) rw? sinh (v2) ” 
+} (1—#?) - i 
sinh (VX(C— )) c+e, 2 sinh (V4) d+d 1, 2 
ae sinh (V0) ( a 3) * San (va ( z tye *) 
ats g 


rd??? 
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sinh (Vr(e- 0) ef od sinh (vXe) =< 2 Lg: 
eS sinh (vr0) ( | ) 2 sinh (v2) (+ . A | ») 
1 
“hk 


2 
2 
(1+#’) — 5. 
3.2. Solution Method 2. The solution is according to the generalized differentiability. Con- 
sider a—level sets of the boundary value problem (3.1), that is 


uy 


y (t) = dy (t) + [P]",¥0) = [6" 9 = hl", (3.12) 


where 
(i? -—1+0,?+1-a], 


a= [es Ea (2)a 
a= [a (258) 03 (®H) 


Also, (i,j) solution means that y is i-differentiable and y’ is j-differentiable (i,j=1,2). 
Using the generalized differentiability and fuzzy arithmetic, for the solutions (1,1) and 
(2, 2), 


pare 
ga 
Q 

l| 


y= Ay +t?—1t+a 
y, (0) =e+ (=) a (3.13) 
y, @Q=dt+ (54) av 

a 


a (3.14) 


(3.15) 


must be solved. 
(1) The solutions (1,1) and (2,2) 


From the solutions of the differential equations in (3.13) and (3.14), the lower and the upper 
solutions of the boundary value problem (3.12) are obtained as 


Vat 


_VJnxt 15 2 1 
t | 1 
+ Coe . X ( a), 


y, (t) = cre 
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=. RO <r eee 2 1 
Yq, (t) = tieY™ + Ge rs w +5 (a —1). 


Using the boundary conditions, c,, Cy, €1, €2 are solved as 


(a (S2)a t= 4 + (C +a-1)) —e-vM (c+ (5) te — <(1 a) 
= evrt — e— VE : 

(e+ (S2)a+ 2-F(1-a)) - (d+ (Sf)a+ 344 (C+0-1)) 
OFS evil _ e—Vit , 

(a (54) at+a+z(C41- a)) see 0 (< (5) a+ —-Fz(a 1) 
aS evr — e- VE 
= eune (z (5) at ~- x(a Je ales at ot 5(P4+1- a)) | 


From this, for the solutions (1,1) and (2,2) the solution of (3.12) is 


wl = [y, 7.) (3.16) 
sinh (VA (¢— 2) oon a 
y(t) = Sah (V3) (c+( ; Jor wa a)) (3.17) 
sinh (V4) d—d ae 
"aan (V30) (a ( 2 Jers es +a-1)) 
ae cae emt 
- sinh (VA(E-1)) 7 (pe ey 
Ya(t) = sin (V%0) (< ( 5 Jo. aa. 1) (3.18) 
sinh (V4) a d—d Dis 
en) (a ( 2 Jers a +1) 
1g AO 
xt yz +x (0 1). 


Proposition 3.1 The solution (3.11) according to the solution method 1 is the same as the 


solution (3.16) according to the solution method 2. 


Proof If a — cut of the solution (3.11) is taken, we have 


y)=y0+ (Po) a, a. =90- (PX 
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That is, the proof is complete. 


Theorem 3.1 The (1,1) solution of the problem (3.12) is a valid a— level set for te [0,2 


satisfying the inequality 


tanh (vat) = (v%¢) ~ (ery) > 0, (3.19) 
sinh (v2) 


The (2,2) solution of the problem (3.12) is a valid a— level set for te [0,4] satisfying the 


tanh (vat) a (v%") ~ (5) <0. (3.20) 
sinh (v2) 


inequality 


Proof If 
S Ya (t) 7 


, 


0 = 
Ba 0, ne <0, 9, () S Fat), y, () 


the (1,1) solution of the problem (3.12) is a valid a— level set. 


If 
a .S ‘ a uw 
Yul) | 9 Walt) 0,4, (t) < Fat), alt) <y, (t) andy” (t) <7, (t), 


0 
Ja "0a =a 
the (2,2) solution of the problem (3.12) is a valid a— level set. 


For the (1, 1) solution, 


dy, (t) sinh (vr (e- t)) 


Ja sinh (vr2) 


(580). 
2 A sinh (Ve 


sinh (Vr(e-8) J 
Ja sinh (v22) ( 


sinh (v22) 
sinh (v8) (a oe 2) ) in 
sinh (v22) 


J.-y,O = C 


pee) (. ” " 
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Also, derivating of (3.17) and (3.18), we have 


is : VX cosh (vx (€— t)) a 2 
Yq (t) eis, (t) = (1-a) oa (Vx) (< c+ ;) 
VX cosh (vxe) 7 2 
ann (V0) (-a+5) 
Then, if 
cosh (vx) (a—a+ 5) > cosh (VrA(e- 0) (e-c+ 5) ‘ 
we have 


From this, making the necessary operations, it must be 


c-¢e+y 


sinh V0) 


Also, again derivating of (3.17) and (3.18), we have 


or sin Vr — 
moa Or S-6 o(° ne) € By :) 


cosh (v22) _ (=) 


tanh (vat) > 


sinh (v2) A 
A sinh (vxt) i 2 
| sinh (v2) (a a 3) oe 


Consequently, the (1,1) solution of the problem (3.12) is a valid a— level set for te [0,2 
satisfying the inequality (3.19). For the (2,2) solution, the proof is similar. 


Theorem 3.2 For any t € [0,¢], the solutions (1,1) and (2,2) of the problem (3.12) are 
symmetric triangle fuzzy numbers. 


Proof Since 


sinh (VA(C-2)) (ete 2 
sinh (v2) ( 2 oo) 


y, 


_sinh (vt) (i44, 2 2) 1 2 
“sinh (VX2) ae ae r d?2 


a on (t) ’ 
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inh (V/A (€ — e- 

way = om (MHD (9 
sinh (vxt) ae ~ th ik 
ee (vX2) ( 2 . 7 r 


= Ya (t) —N (t) ) 


the solutions (1,1) and (2,2) of the problem (3.12) are symmetric triangle fuzzy numbers for 


any t € (0, 4). 
(2) The solutions (1,2) and (2,1) 


Using the generalized differentiability and fuzzy arithmetic, the solution of (3.15) is 


iv Ol = |v, © 5.0], (3.21) 
1 2 1 
y,, (t) = ce + ce" — cg ain (vat) — €4 COS (vat) il ty (1-a), 
1 2 1 
Yq (t) = ev ee c3 sin (vat) + €4 COS (vat) x 2h (l—a). 


From the boundary conditions, the coefficients c1, co, c3 and c4 are obtained as 


1/2, 2 , dtd —Vrxe (2 | tte 

(40 et ee) e wt 

C= ’ 
eV — e-VE 


(24 He) (42 2 + 4) 


2 EVM ene ) 
(1 — a) @ = cE #4) nT 
eo in(Vae) WX 


1 G-e 


asa 


Theorem 3.3. The (1,2) solution of the problem (3.12) is a valid a— level set for te (0, ¢] 


satisfying the inequalities 


eg ge C—C 
wm) sin (VAe) + (5 * -*) cos (VXt) — : = 0, (3.22) 
on cos (vat) - (5 + aa sin (vx) > 0. (3.23) 


The (2,1) solution of the problem (3.12) is a valid a— level set for te (0, ¢] satisfying the 
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inequalities 
24 te, ded 
anaes ma eS | (vre) + (5 + ) cos (vat) SES Gy (3.24) 
sin (v22) r 2 r 
2, te, dd a 2 
peal were ag (vat) = (5 + -*) sin (vx) <0. (3.25) 
sin (v2) r 2 
Proof If 
dy , (t) dy, (t) ; ; " if 
=o = <a <y 7 ‘ 
Fat aq <8 Ya) So 4), Yq) S Go (4) and F(t) Sy, (4), 


the (1,2) solution of the problem (3.12) is a valid a— level set. 
If 


dy , (t) dy, (t) : ; ” i 
=o o <qy y < y < 
Fa Bq <8 Ya) So (4), Ga (t) Sy, (4) and Fa (t) Sy, (4), 


the (2,1) solution of the problem (3.12) is a valid a— level set. 
For the (1,2) solution, from 


dy, (¢) 
Oa 


Oa (t) 


0 
abe, Oa 


<0 and y, (t) —y, (t) = 0, 


it must satisfies the inequality (3.22), from 7, (t) — vy, (t) > 0, it must satisfies the inequality 
(3.23). Also, from y(t) — 7. (t) > 0, it must be 


12) avi) + (F458) n (v0 


Then, the (1,2) solution of the problem (3.12) is a valid a— level set for t€ [0, 4] satisfying 
the inequalities (3.22) and (3.23). For the (2,1) solution, the proof is similar. 


Theorem 3.4 For any t © [0,], the solutions (1,2) and (2,1) of the problem (3.12) are 


symmetric triangle fuzzy numbers. 


Proof Since 


1 2 
y, t= ce! + cge-V** — 7? - — = y(t), 


y, © —y, (0) = egsin (VXt) + ex 008 (Vat) — : (1—a) =%, (t)- 7, (t), 


the solutions (1,2) and (2,1) of the problem (3.12) are symmetric triangle fuzzy numbers for 
any t € (0, 4). 
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Example 3.1 Consider the fuzzy boundary value problem 


uw 


y (t)=y(t) + [P]" 9) = [0]*,¥Q) = [1I", (3.26) 


and [0]* = [-1+a,1-— al, [1]*° = [a, 2 — a]. Then, the (1,1) and (2, 2) solutions of the problem 
(3.26) are 
_ 2asinh(1—t) | (2a +2) sinh (¢) 


A Cs ar TT Tc aa (3.27) 

= _ (4—2a)sinh(1—t) | (6—2a)sinh(t) 5 
Ja (t) = sinh (1) "sinh (1) ere 2) 
iy Ol = |v, Fa ()]- (3.29) 


The (1,2) and (2,1) solutions of the problem (3.26) are 


i) = (—S)¢ ! (-S)q (A) sine (3.30) 


—2(1—.a)cos(t) -—#-—a-1, 


Ya (t) = (>) e4 (= ) e*4 () sin (t) (3.31) 
+2(1—a)cos(t) -—? +a—3. 


iy Ol = |v, Fa ()]- (3.32) 


The (1,1) solution is a valid a— level set for te [0,1] satisfying the inequality tanh (t) — 


(cau) > 0, the (2,2) solution is a valid a— level set for t€ [0,1] satisfying the inequality 


tanh (t) — (S44) <0. 


Figure 1. Graphic of the function tanh (t) — (a) 
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According to Figure 1, the (1,1) solution is a valid a— level set for t > 0.5 and the (2, 2) 


solution is a valid a— level set for t < 0.5 

The (1,2) solution is a valid a— level set for t€ [0,1] satisfying the inequalities ee + 
2cos(t)-1> 0 and Sa —2sin(t) > 0. The (1, 2) solution is a valid a— level set for té [0, 1] 
satisfying the inequalities ma sin (t) + 2cos(t) —1> 0 and mD cos (t) — 2sin (t) < 0. 


Figure 2. Graphic of the function ee + 2cos(t)—1 


Figure 3. Graphic of the function a — 2sin (t) 


According to Figure 2 and Figure 3, the (1,2) solution is a valid a— level set and the (2, 1) 
solution is not a valid a— level set. 

Also, the solutions (3.27)-(3.29) and (3.30)-(3.32) are symmetric triangular fuzzy numbers 
for any t € [0,1]. 
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Figure 4. Graphic of (3.27)-(3.29) for a = 0.5 
Blue + J, (t), Red—y (t), Green + y, (t) = y, (t) 


w 


Figure 5. Graphic of (3.30)-(3.32) for a = 0.5 
Blue + y, (t), Red —y, (t), Green > 9; (t) = y, () 


§4. Conclusion 


In this paper, a problem with fuzzy forcing function and fuzzy boundary values is investigated. 
The problem is solved by two different solution methods. It is found that the solution of the 
problem according to the solution method 1 is the same as the solution according to the solution 
method 2 for the solutions (1,1) and (2,2). It is shown whether the solutions (1,1), (2,2), (1,2) 
and (2,1) are valid a—level sets or not. Also, all the solutions are symmetric triangle fuzzy 
numbers for any t € [0,4]. Example is solved. It is shown whether the solutions are valid fuzzy 
functions or not. Graphics of solutions are drawn. The solution method 2, that is, solution 
method according to the generalized differentiability is more useful than the solution method 
1. Because, it is found the wider solutions of the problem with the solution method 2. 
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The Variation of Electric Field 
With Respect to Darboux Triad in Euclidean 3-Space 
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Abstract: In this paper three electric fields are described via Darboux triad components 
in Euclidean 3-space. Later variations of three cases of electric field with respect to Darboux 
triad are studied. Finally Lorentz force equations are presented via electromagnetic magnetic 


curves with respect to Darboux triad in Euclidean 3-space. 
Key Words: Geometric phase, Darboux frame, electric field. 
AMS(2010): 53A35, 53B30, 78A05. 


§1. Introduction 


The geometric phase is described as the angle of rotation a light wave travelling in optic. The 
phenomenon of a geometric phase have many applications in condensed-matter physics, optics, 
particle physics, gravity, cosmology, chemical physics and mathematics [1-6]. The geometric 
phase is connected with parallel transport of the polarization along curved light [7-9]. 


Berry studied adiabatic phase and Pancharatnam’s phase for polarized light [10]. Recently 
numerous authors presented the the electric field variation of along an optical fiber [11-14]. 


Balakrishnan et al. presented anholonomy density via Frenet triad in Euclidean 3-space 
€3 [15]. Three geometric phases and parallel transports for numerous frames have been investi- 
gated by Giirbiiz in [16-20]. Balakrishnan introduced geometric phase for first class associated 
with some solitons for Darboux triad in €° [21]. New classes associated with the nonlinear 
Schrodinger NLS equation for Darboux triad in €° have been given in [22]. 


The electric polarization theory contains the geometric phase phenomenon [23]. Mukunda 
and Simon showed that the unit electric vector field E is written via the principal normal vector 
field N and the binormal vector field B of the Frenet triad {T, N, B} in Euclidean 3-space [24]. 
In this paper we express three electric fields via Darboux triad apparatus. Later evolutions of 
three electric fields are studied via Darboux triad in €°. Eventually Lorentz force equations are 


obtained via electromagnetic curves with respect to Darboux triad in €%. 
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No.202019016). 
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§2. Preliminaries 


Let T', be a curve on a connected surface S with the arc length o in €3. Apart from Frenet 
triad, at every point of curve, there is a Darboux triad {t,g,n}. t is the tangent vector, n 
is the normal of surface and g = t x n. The spatial evolution of the Darboux triad {t,g,n} is 


given by [25] 
to 0 Ke Ko t 
Zo = — 9) 0 7) g (1) 
Ne — KS) —7{9) 0 


KO) is the geodesic curvature, the normal curvature is K) and 7) is the geodesic torsion of 


the curve I. The time evolution of the Darboux triad {t, g,n} is given by 


£3} 0 K6o) (2) t 
g. |=| —n? o x” g (2) 
Nn, 16°) —7{0) 0 n 


where u denotes time and t,, = ot 


A magnetic field is a closed 2-form F in €?. The Lorentz force © of a magnetic background 
(E€%,(,)) is a (1,1) type skew-symmetric tensor and it is described as 


F(z, y) a (z, y) 


x,y € x(E%). A smooth curve [ in (€3, (,)) is described as a magnetic curve of the dynamical 
system connected with the magnetic field F if its velocity vector field satisfies the following 
differential equation I.,, = ®([I,). Divergence free vector fields and magnetic fields are one to 
one correspondence, the Lorentz force ® concerned with the magnetic field M [26], [27] 


O(r) =MaAz. 


§3. Geometric Phase for First Case of Electric Field with Darboux Triad in €° 


Balakrishan introduced first frame {P,,P2,P3} and first transformation € of curve evolution 
concerned with the NLS equation with respect to Darboux triad in €° as following [21] : 


g +in ;¢¢ do’ * g- im -if? Ido" 
P, = t,Po= bl tte ps e a (3) 
J2 J2 
fn, EEO ee. (4) 


V2 


The spatial evolution of the first frame {P1, Po, P5} is given by 


Pig =€*P2+€P5, Poo = —€Pi1, P35, = —*Pi (5) 
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where €* is the conjugate of €. Also temporal evolution of {P;, P2, P3} is 


Pi, ty = —\"P2 — \Pt (6) 


I 


Py. AP, + iZP2 (7) 


I 


where Z(o, wu) is a real function. From Poy, = Pau, it can be obtained: 
Ty = 1rAE* — id*E. (8) 


where 
AD, dodu = (780) — Ths) )dodu 


is first anholonomy density measure for polarization plane of linearized light wave travelling 
along optic fiber in €3 [21]. 
ae I Baan (9) 
V2 


satisfies Eqs.(6), (7) and (8). The time evolution of the Darboux triad is given by 


t, = 6?) xt=rgt+uwn (10) 
Su = f) xga-rt+ 7n (11) 
n, = &) xn=—wt— Tog (12) 


where of?) =( (+4 Big+Cin), r= C\,w = —B,. Using Eq.(4) and Eq.(9), Z> = KO r— Kw, 
The time evolution of Darboux triad for first class can be written by 


ty =rg+wn (13) 
1 
2, =—rtt+ (f tS) do —T)n (14) 
O1 ? 
n, = —wt — (f ido —L)g (15) 


and anholonomy density 
AD, (o,u) = -I, = —rrK9) + wr) 


for first class. Total phase P for first class with respect to Darboux triad in Euclidean 3-space 


-{ i, Tyilodu= | / (t,t, X t,) dodu 
U1 oo U1 90 


U2 O1 
= i if (—rK6s) + wk) )\dodu 
Uy oo 


is given by 


Nu 
I 
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Also [22] 
Pir = —i€%P2+i€P5 
Po, = —i&Pi+TZP2, 
PS, = ws —TP3, T =iée". 


From Piys = Pigy and Pous = Pasy, the NES equation system 


l 


ei iba +i\é|?& 
& itso —il€|? é. 


I 


is obtained. 


A optical fiber can be described by the curve [';(c) on any surface with respect to Darboux 


triad in €3. The change of the electric field E, can be written by 
Eio = pit + yog + Ysn. 


Case 1. Assume that 
(E,,t) = 0. 


Using Eq.(16) and Eq.(17), it can be obtained 
1 = —Kg (E1,g) — kn (Ei,n) . 
When no various loss mechanism along the optic fiber, 
(E,, E,) = const. 
Using Eq.(16) and taking derivative with respect to o of Eq.(19), it can be derived 
y2 (Ei, 8) = —93 (Ei, n). 
Via Eq.(20), it can be obtained 


y2 = @(Ei,n), 93 = —(Ey,g) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


The evolution for the polarization of light wave travelling from the point ['(a9) to the point 


T1(o1) along the Ty = T4(o) curve with respect to Darboux triad is given by the evolution of 


the electric field E,. 


Consider (E,,g) # 0, (Ei,n) 4 0. Substituting Eqs.(18) and (21) in Eq.(16), the change 


of the electric field E; is written by 


Ey, = (—Kg (E1, 8) — Kn (Ei, n))t+o@ (E1,n) g—w (Ei,g)n 


(22) 
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where @ is a parameter. Using Eq.(20) for @ = 0, Eq.(22) is rewritten by 
Eto = (—Ky (E1, 8) — Kn (Ei, n))t (23) 


The Fermi-Walker derivative of the electric field E, with respect to Darboux triad in €? is 
given by 
seal Dee = E,. — (t, E;) te + (to, E;) t. (24) 


The electric field E; is the Fermi-Walker parallel transport if and only if 
DFWE,, = 0. (25) 

Using Eqs.(17), (24) and (25) it can be obtained 
Ei, = (to, E;) n. (26) 
The electric field vector E, with aid of the Darboux triad apparatus g and n is expressed 

by 

(g + in) g—in 
v2 V2 - 


where E, Ei = 1 and |Q(c)|? + |Q*(c)|? = 1, E* is complex conjugate of Ey. 


O1 
p= | 1)do 


is the change phase of the polarization light injected into this fiber with respect to Darboux 
triad in €°. 


Eyo(o) = Q(a) +0*(c) 


(27) 


Q(o) = eI 84" O(a) 


OF (a) = Pee he 7§do" O* (gp) 


with the polarization coefficients are 


Also via Po, P§, Q(00) and Q* (a0), the electric field E,(c) is expressed as 
E, (0) = P2Q(00) + P3O*(00) (28) 


Respectively, taking derivative with respect to o and the time u of Eq.(28), the spatial and 
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temporal evolutions of the electric field E; for Darboux triad are derived as following: 


Ey, => Po,Q(c0) + P32" (0) 
Fux = Poy2(d0) a P52" (00). 


From compatibility condition Ey, = Ej, ,, the nonlinear Schrodinger NZIS equation 
system. 


The Lorentz force equation 6“ of the electric field vector Ey is given by 


OOF, =F, =M” x By (29) 

and 
(@E,,t) 3s (Ey, 8), (®E,,g) me (E:, 8g) (30) 
(@E,,n) ss (E:, on). (31) 


The trajectory of travelling particle along the magnetic field M“ with respect to Darboux 
triad in €? is described as electromagnetic trajectory. If DEM” curve follows the magnetic 
trajectory, it is described as the Darboux electromagnetic curve in €°. With the help of Eqs. 
(30), (31) the Lorentz force ®* equations in the Darboux force equations of the DEM™ curve 
of the [; are given by 


6) (t) 0 — 4) — 4) t 
SH (g) | = Ke 0 -—w g (32) 
(4) (n) Ko ow 0 


DEM”) curve of the [; is a magnetic trajectory of the magnetic field M“ divergence free 
field iff M) is given by in the following 


M® = —ot + Kn@—KgN 


§4. Geometric Phase for Second Case of Electric Field with Darboux Triad in €? 


Respectively, the second frame {Q1, Q2,Q3} and second transformation ¢ associated with the 
NLS equation via Darboux triad is given by [22] 


Qi = &, (33) 


Q, t+ Un gy KS) do’ _t—-m 


* if KO do’ 
pre ee (34) 


(s) 


d = (—Kg 


+ iTg’) ise nao’ 
e€ 


a (35) 
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Using Eqs.(33) and (34) the spatial evolution of the frame {Q1, Q2, Q3} is given by 


Qic = $'Q2.+9Q3 


Qo, = —¢Q, 
Q3, = -¢'Qi 
st, (—K) -i7) -is? ldo 
where 0) = Safa i 
Consider 
Qiu = Su = 42Q2 + b2Q3 + e2Qi (36) 
Qou = h2Q2+ foQ3+VQi. (37) 


From (Qix, Qi) = 0 > cz = 0, (Qiu, Q2) = b2, (Q2u, Qi) = 0 > bo = —-9, 
(Qou, Q2) am fe a 0, (Q3.,, Q2) oa —hg => he = —fs and ag = —v*. 
Eqs.(36) and (37) are rewritten by 


Qiu = Su=—V"Q2 — IQ) (38) 
Qa = VQi+i7Qe2 (39) 
with J(c,u) areal function. From Qayo = Qaoy, the followings are obtained 


bu 
Io 


l| 


—U, +iJSo 


I 


Wg" — i*¢. (40) 


O1 ; 
When t and n rotates around g with KS) (a), a geometric phase P = f[ KS) (a)do arises 


20 
between t, n and corresponding nonrotating Darboux triad in €°. 


When the linearized light wave travelling moves from wu, to ug along the curve in optic fiber, 
u2 
a geometric phase P = [ KO?) (u)du arises between natural Darboux triad and nonrotating 


U1 
Darboux triad in Euclidean 3-space. The rotation angles of polarization plane can be given by 


Pi = KI (a,u)Ao + 6) (o + Ao, u)Au 
Py = 6) (o,u)Aut KO (o,u+ Au)Ao 


The phase difference are given by 6P = P; — Pz = ADa(o,u)AcAu. 
AD»2 = (xi) - KO) is second anholonomy density measure for polarization plane of linearized 
light wave travelling along optic fiber for second case in €°. Also 

(J a5 iw) fc Kido 


v= re a (41) 


satisfies Eqs.(39) and (40). The time evolution of Darboux triad for second class is given by 
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[22] 
ty = (?xt=—igtKOn (42) 
g. = (?xg=lit+wn (43) 
n, = of?) xn = —Kt — wg (44) 


where of?) = Ast — Ko) p + Con, | = —Co, w = — Ao. 
Using Eqs.(35), (40) and (41) it can be obtained 


Jo = (TOL + Kw). (45) 


From Eqs.(34), (39), (42), (43) and (44), the time evolution of Darboux triad for second 
class with Eq.(43) is given by 


o1 
tule +(f KS)do — J)n 


O1 : 
ist / do’ — T)t— we 


and the anholonomy density AD2(o,u) = —J, = (71 + KS w) for second class. Total phase 
P for second class with respect to Darboux triad in €° is given by 


U2 O11 U2 o1 
P= -| / Jadodu = | / (791+ 69) w)dodu 
U1 oo U1 00 
U2 O1 
[of (eee x eu) dow, 
U1 (exe) 


The quantum geometric phase for second class of curve evolution with respect to Darboux 


triad in €° is obtained Bt. 2a 3es 
Pail dos. (Qou,Q%) du. 


l| 


Also [22] 


Qiu 
Qo 


I 


—16,Q2 + 196Q3, Qou = —i6cQi + FQz, 
16,Q1 — JQ3, JT =igoe" 


l| 


From P1uc = Piou and dbouc = db2cu; the NLS equation 
bu =itboo til oo 


is obtained. 

A optical fiber can be described by a curve ['2(c) with respect to Darboux triad in €°. The 
direction of electric field Ez is given by the direction of the state of the linearly polarized light 
wave injected to the fiber with respect to Darboux triad in €°. The change of the electric field 
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E»2 with respect to Darboux frame in €° can be given by 
Eo, = Cit + Cog + Can. (46) 


Case 2. Assume that 
(E2,g) =0. (47) 


Using Eqs. (46) and (47), it can be written by 
(2 = —Kig (Eg, t) — Ty (E, n) (48) 


Consider 
(Eg, E2) = const. (49) 


Taking derivative with respect to o of Eq.(49), the followings are obtained 
G1 (Ba, t) = —C3 (Ba, n) (50) 
G1 = x (E2,n), ¢3 =H (Ez, t) (51) 
where y is a parameter. 


Using Eq.(20) and (E2,t) 4 0, (Eo,n) 4 0. Substituting Eqs. (48) and (51) in (46), the 
evolution of the electric field vector Ez with respect to Darboux triad is given by 


Ex, = x (Ez, n) t+(—kg (Eo, t) — Ty (E2,n))g — x (Ee, t)n (52) 


Via Eq.(52) for x = 0, 
Ey, = (—Ky (Ep, t) — Tg (E2,n))g (53) 


The modified Fermi-Walker derivative for the electric field vector Ez with respect to Dar- 
boux triad for second class is described by 


Dm PW Boo = Exo — (8, E2) 8c + (Bo, Ez) & (54) 
The electric field Ez is the modified Fermi-Walker parallel if and only if 
DmFWE,, = 0. (55) 
Via Eqs.(47), (54) and (55), one obtains Ey, = (g,, Eo) n. 
The electric field vector Eg with aid of the Darboux triad apparatus t and n can be 


expressed by 
(t + in) t—in 


Qo Ce 


Eo(c) = T(c) (56) 


where E2E3 = 1 and IT (o)|? + |Y*(0)|? =, 
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Here T(c) and Y*(c) are 
—is71 nb) do’ 


L(o) = eS KP de V(g), T*(o) =e T*(o0) (57) 


and the polarization coefficients are 


T(a0) = (SE) Baloo), T(e0) = (*) Baloo) (58) 
Via Eqs.(34) and (57), Eq.(56) is re-expressed by 
E2(7) = Q2T (a0) + Q3T* (a0). (59) 


Respectively, the spatial and temporal evolutions of the electric field Ez for Darboux triad 


are derived as following: 


Ey, = Qo>T (a0) F Q3,T" (70) 
EB. = QouT (90) + Q3,,T* (90). 


From compatibility condition E2,, = E2,,, the NES equation system connected with the 
electric field Ey is derived. 


Geometric phase for polarized light injected into a fiber with respect to Darboux triad for 


o1 
P =) KSdo . 


Consider the Lorentz force equation (9) for second case of the electric field vector 


second case in €° is given by 


& Fy = Ey, = M9) x Ep (60) 

and 
(@E»,t) a age (Es, 8) (®E,,g) = (B;, &g) ; (61) 
(®B,,n) =a (Eo, 8n). (62) 


The trajectory of travelling particle along the magnetic field MY with respect to Darboux 
triad is described as the electromagnetic trajectory. If DEM) curve follows the magnetic 
trajectory, it is described as the Darboux electromagnetic curve. With the help of Eqs. (61) 
and(62), the Lorentz force ®9 in the Darboux triad of the DEM"? curve of P'2 are given by 


#9) (t) 0 —«{9) —-x t 
&9(g) | = Ki) 0 ine g (63) 
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Via Eq. (63), the vector field divergence free M‘) is given by 


M = yg + Tgt—KgN. 


§5. Geometric Phase for Third Case of Electric Field with Darboux Triad 


The third frame {R,, Ro, R5} and the third transformation w for third class of curve evolution 


concerned with the NLS equation with respect to Darboux triad in €° are given by [22] 


R, = Qn, (64) 
t+ig 5 fer pSdo’ t—i £71 nde’ 
Ry = TBI Pde’ Ry a Bi Hae (65) 


ae) 


x 


(Ki? + iri?) $7 nhs) do” 


= 66 
vy we (66) 
Using Eqs. (65) and (66), the spatial evolution of {R,, Re, R43} is given by [22] 
Rig = —*Re — PR3, Reo = ~Ri, Ro, = Y*Ri (67) 
sie (665) —ir{)) —is? wh do’ 
where w* = ary ame F 
Consider 
Ri, =n, = a3Re + bsR5 + c3R1, (68) 
Ro, = h3Ro + fsR5 + 7R4. (69) 
From (Riu, R1) =0=> Cc = 0, (Riu, R2) = bs, (Rou, R1) =n=> b3 =-Y”, 
(Rau, Re) = fs = 0, (Riu, R35) =a3> n* = a3, (R3,,, Re) = —hg > hg = — fz. 
Eqs. (68) and (69) can be rewritten by 
Riu = Ny = —7*Re a nV5, (70) 
Ro, = nR, + iLRe (71) 
with £L(c,u) is a real function. From Rous = Raou the followings can be derived by 
Wu =o + Lo, (72) 
Ly = int — inu™. (73) 
1 
When t and g rotates around n with Ke) (o), a geometric phase P = [ Ko do arises 


20 
between t, g and corresponding nonrotating Darboux triad in €?. When the linearized light 


wave travelling moves from u; to ug along the curve in optic fiber, a geometric phase P = 


U2 
a Ko? du develops between natural Darboux triad and nonrotating Darboux triad in €°. The 
U1 
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rotation angles of polarization plane can be given by 
P= Ki) (a,u)Ao + Koo + Ao, u)Au 


P2 = KO (a, u)Au+ Kis) (g,u+ Au)Ao. 


Phase difference are given as /P = P, — Pz = AD3(c,u)AoAu, where AD3 = (xs) - KO) 
is third anholonomy density measure for polarization plane of linearized light wave travelling 
along optic fiber for third class in €3. Also 


© + iz) i fr (gq : 
——_“ _—“e’ Reg G0. (74) 
V2 


satisfies Eqs.(70), (71) and (73). The time evolution of Darboux triad is given by [22] 


ty = 65" x t Kf) g — jn, (75) 
Lu = f xg =K Ot — zn, (76) 
ny, = of) xn=jt+ zg (77) 


where of) = A3t+ Bsgt+ Kn, z = —As3, j = Bs. Using Eqs. (66), (74) it can be obtained 
Le jr - KhS) z, (78) 


The time evolution of Darboux triad for third class connected with the NLS equation is 


given by 
n, = jt+ zg (79) 
o1 
ee ee (f KQdo' — Lg (80) 
O1 ‘ 
gy = (£- KS) do )t — zn. (81) 


The anholonomy density AD3 for third class with respect to Darboux frame in Euclidean 
3-space: 
AD3(0,u) = Ly = Kz — ari. (82) 


and the total phase P for third class with respect to Darboux triad in €° is given by 
U2 O1 
-{ ) L,dodu 
U1 (exe) 
U2 O1 U2 O1 
= / i (Kz — jt) )dodu = if / (n, ng X n,) dodu. 
UL 00 Uy [oxe) 


The quantum geometric phase is given 


aa ra) u2 . 
Pail dos, [ (Ry, Ri) du. 


P 


I 
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A optical fiber can be described by a curve ['3(o) with respect to Darboux frame in €°. 
The direction of electric field E3 denotes the direction of the state of the linearly polarized light 
wave injected to fiber with respect to Darboux frame in €°. The change of the electric field E3 
with respect to Darboux triad can be written by 


E3o =7™t+ meg + 73n. (83) 
Case 3. Assume that 


(BE, n) = 0. (84) 


From Eq.(84) 
(E3,, n) = (E3,knt+T,g) 


13 = Kn (E3,t) + Ty (Es, g) (85) 


Also 
(E3, Es) = const. (86) 


Using Eq.(83) and taking derivative with respect to o of Eq.(86), it can be obtained 
TY (Es, t) = —T2 (Es, g) (87) 


m™ = €(E3,g), 12 = —€ (Es, t) (88) 


where € is a parameter. The evolution in the polarization of light wave travelling from the point 
T'3(09) to I'3(01) along curve with respect to Darboux triad is given by the evolution of the 
electric field E3. (E3,t) 40, (E3,n) # 0. Substituting Eqs. (85) and (88) in (83), the Eq.(83) 
is rewritten by 

Ego = € (Es, 8) t — (Ba, t) 8+(tin (Es, t) + 7y (Bs.g))n (89) 


Via Eq.(89) for « = 0, 
E3, = (Kn (E3,t) + 7, (E3,g))n (90) 


The modified Fermi-Walker derivative for the electric field E3 with respect to Darboux 
triad for third class is described by 


2mEW Eso = Ezq — (n, Ez) ng + (i, Ez) n (91) 
The electric field Ez is the Fermi-Walker parallel if and only if 
DmFWR,, = 0. (92) 


Via (84), (91), (92), one obtains E3, = (ng, Ez) n. 
The electric field vector E3 with respect to the Darboux triad apparatus t and g can be 

written by 

t —ig 


(t + ig) o 
+ U*(o) : 
V2 


E3(c) = H(c) 5: | 


(93) 
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where E3E% = 1 and |X(c)|? + |E*(o)|? = 1. Here 


~ , 
-isti nde 


E(o) = ef 40" y(g9), U*(a) =e E*(ao). (94) 
The polarization coefficients are 
t+ig\* 
U(oo) = E3(a 
(00) = (S738) Baloo) 
t— ig e 
b* (09) = E3(d0). 
(00) = (SE) Baloo) 
Eq.(93) is re-expressed as the following 
E3(o) = R2U(o0) + R3=*(o0) (95) 


When taking derivative with respect to o and the time u of Eq. (95), the spatial and 
temporal evolutions of the electric field E3 for Darboux triad are derived as follows 


E35 = ReoX(o0) + R5,%* (00) 
E3y, Roud(o0) + R35," (a0) 


From compatibility condition E3,,, = E3,,, the nonlinear Schrédinger equation NLS 
system connected with the electric field E3 is obtained. 


Bei KS do" 


is the change phase of the polarization light injected into a fiber for third case of the electric 
field with respect to Darboux frame in €°. Consider the Lorentz force equation ®) for third 
case of the electric field vector 


@%E, = Es, =M™” x Es, 
(@™Es,t) = — (Bs, Ot), (OEs,g) = — (Es, Og) 
(@ Es, n) = — (E53, 8 n), 


The trajectory of travelling particle along the magnetic field M“” with respect to Darboux 
frame is described as the electromagnetic trajectory. If the curve DEM follows the magnetic 
trajectory, it is described as the Darboux electromagnetic curve. With the help of Eq. (61) 
the Darboux Lorentz force equations along the optic fiber for third case the electric field are 


obtained 
#(") (t) 0 —e Kh) t 
(") (g) = € 7) 0 g (96) 
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DEM" curve of the Ts is the magnetic trajectory of the magnetic field M(”) iff the vector 
field divergence free M“ is given by 


M(”) = —«{S) g — entr)t, 
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Abstract: Recently, a novel invariant is considered, which is the Nirmala index defined as 
the sum of the square root of the degrees of the pairs of adjacent vertices. In this paper, 
we introduce the first and second inverse Nirmala indices of a graph and compute exact 
formulas for certain nanostructures. 
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§1. Introduction 


Let G be a simple, finite, connected graph with the vertex set V(G) and edge set E(G). The 
degree dg(u) of a vertex u is the number of vertices adjacent to u. The additional definitions 
and notations, the reader may refer to [1]. 

A molocular graph is a graph in which the vertices correspond to the atoms and the edges 
to the bonds of a molecule. A topological index is a numeric quantity from structural graph 
of a molecule. Several topological indices have been considered in Theoretical Chemistry, and 
have found some applications, especially in QSPR/QSAR study, see [2, 3, 4]. 

In chemical science, numerous vertex degree based topological indices or graph indices have 
been introduced and extensively studied in [4, 5]. 


The Sombor index was defined by Gutman in [6] as 


SO(Q)= D2 Valu)? + da(v)?. 


uvE E(G) 


Recently, some Sombor indices were studied in [7, 8,9, 10, 11, 12, 13, 14]. 
In [15], Kulli introduced the Nirmala index of a graph G and it is defined as 


NG@= Yo Vdelu)+de(v). 


uv€ E(G) 
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We now define the first and second inverse Nirmala indices of a graph G as 


IN(G) 


I 


oe laa J wai] : 


1 i is 
IN,(G) = > laces 7 ao 


uve E( 


In this study, we compute the first and second inverse Nirmala indices for four families of 
dendrimers. For dendrimers, see [16]. 


§2. Results for Porphyrin Dendrimer D,,P,, 


We consider the family of porphyrin dendrimers. This family of dendrimers is denoted by D,, Py. 
The molecular graph of D,,P,, is shown in Figure 1. 


Pp GENS I 
as agp VG Ry 
3 etd, load es 


gist tn 7 i a ae 3 


yt fe <> 
Wily “a 


Figure 1. The molecular graph of D,P,, 


Let G be the molecular graph of D,,P,,. By calculation, we find that G has 96n—10 vertices 
and 105n — 11 edges. In D,, Py, there are six types of edges based on degrees of end vertices of 
each edge as given in Table 1. 


dg(u),de(v) \ w € E(G) | (1,3) | 4) | (2,2) (2,3) | (3,3) | 3,4) 
Number of edges 2n 24n | 10n—5 | 48n—-—6 | 18n 8n 


Table 1: Edge partition of D,,P,, 


In the following theorem, we compute the first and second inverse Nirmala indices of D,, Py. 


Theorem 2.1 Let D,P,, be the family of porphyrin dendrimers. Then 


4 2 4 
INi(DnPn) = (+ 12V5+ 10448? + 3¥ + en 5 gv 


V3 v6 V3 V3 v6" 


_ 48 V6 v3 , 16Vv3 V6 
IN2(DnPo) = (v34+ S104 48%2 +13 )n 5 6 
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Proof From the definitions and by using Table 1, we deduce 


3 1 


nf 1 2. 
1 172 1 173 ta 

+ (48n — 6 + 13n] = 4 +8 
a 5+5| nl; | nle+ 


V5 J/2 Nn — 


1 172 i. 172 2 
INi(Dn Pn) =2n|> +5] +24n/> +2) + (10n—5)|54 


4 
= | —+12V54+ 10+ 48 + 13 
(= V6 


and 


INg(DyPa) =2n[7 +5] pain] +3] + (10n 5/5 | 
1177? 11)? 11)? 
+ (48n [5 +5| +13n[5 +5] ran] +3 | 
2 48. v6 v3 , 16V3 v6 
=(vV5+ 10+ +s sn 5 6 


§3. Results for Propyl Ether Imine Dendrimer PETIM 
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We consider the family of propyl ether imine dendrimers. This family of dendrimers is denoted 


by PETIM. The molecular graph of PETIM is depicted in Figure 2. 
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Figure 2. The molecular graph of PETIM 


Let G be the molecular graph of PETIM. By calculation, we find that G has 24 x 2” — 23 
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vertices and 24 x 2” — 24 edges. In PETIM, there are three types of edges based on degrees 
of end vertices of each edge as given in Table 2. 


da(u), da(v) \ We E(G) (1, 2) (2, 2) (2, 3) 
Number of edges 2x2” | 16x 2”—18 | 6x 2”-6 


Table 2: Edge partition of PETIM 


In the following theorem, we compute the first and second inverse Nirmala indices of 
PETIM. 


Theorem 3.1 Let PETIM be the family of propyl ether imine dendrimers. Then 


IN\(PETIM) = (V6+16+V30)2” — (18 + V30), 
IN)(PETIM) = (2 16+ VE) a" (1s we). 


Proof From definitions and by using Table 2, we derive 


mm(periay axel g]: 400.2 wl 3] +6x27-of+ 3) 
= (V6 + 16 + V30)2” — (18 + V30), 

IN,(PETIM) cee) pee Ee 18)|5 | i | exar-o[et]” 
= (22 + 16 4 ee a (18 | ~). 


§4. Results for Poly Ethylene Amide Dendrimer PETAA 


We consider the family of poly ethylene amide amine dendrimers. This family of dendrimers is 
denoted by PETAA. The molecular graph of PET AA is presented in Figure 3. 


red yayt? 


Figure 3. The molecular graph of PETAA 
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Let G be the molecular graph of PETAA. By calculation, we find that G has 44 x 2” — 18 
vertices and 44 x 2” — 19 edges. In PET AA, there are four types of edges based on degrees of 
end vertices of each edge as given in Table 3. 


da(u), da(v) ‘ UWE E(G) (1, 2) (1,3) (2, 2) (2,3) 
Number of edges 4x2” | 4x 2"—2 | 16x 2” —8 | 20x 2”-—9 


Table 3: Edge partition of PETAA 


In the following theorem, we compute the first and second inverse Nirmala indices of 
PETAA. 


Theorem 4.1 Let PETAA be the family of poly ethylene amide amine dendrimers. Then 


IN\(PETAA) = (3 = + 16+ 20)» (= 3+ 8), 


v2 V3 v6 V3 v6 
IN}(PETAA) = (2 2vB+16+ WEE) on (v3+8+ 2). 


Proof By using definitions and Table 3, we obtain 


i 


IN(PETAA) =(4x2")|7 45] + (4x2"—2)[7 45] +(x 2"-9)] 545] 
+ 0x2 —9)/5+5| 


4/3 8 16-+ 2)» (4 we). 


(tat 


T 


and 


IN2(PETAA) = (4x 2") cn le (ase 2)|3 As Al 


1 1717? 
+ 20x 2 —9)/5+5| 


= (42 + 2/3 +16 4 ee 2" (v3+8+ 2). 


§5. Results for Zinc Prophyrin Dendrimer DPZ,, 


We consider the family of zinc prophyrin dendrimers. This family of dendrimers is denoted 


by DPZ,,, where n is the steps of growth in this type of dendrimers. The molecular graph of 
DPZ,, is shown in Figure 4. 
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Figure 4. The molecular graph of DPZ, 


Let G be the molecular graph of DPZ,,. By calculation, we obtain that G has 56 x 2" — 7 
vertices and 64 x 2” — 4 edges. In DPZ,,, there are four types of edges based on degrees of end 


vertices of each edge as given in Table 4. 


da(u), da(v) ‘ UU Ee E(G) (2, 2) (2, 3) (3, 3) (3,4) 
Number of edges 16 x 2”—4 | 40 x 2" —16 | 8x 2"+12 4 


Table 4: Edge partition of DPZ,, 


In the following theorem, we determine the Nirmala index and its exponential of DPZ,. 


Theorem 5.1 Let DPZ,, be the family of zinc prophyrin dendrimers. Then 


(16 40/5 + 2) (4 16V5 12/2 1), 


een v6 " ve VB VW 
runny = (ne oSB)n (oo BE 2) 


Proof From definitions and by using Table 4, we deduce 


wR 


1 
{1/7 1 178 
IN,(DPZp) = (16 x 2” — 4) E + ;| + (40 x 2" — 16) E + ;| 


3.3 3. OCA 


< (oe tr (8 E28) 


i. 472 {. 473 
+ (8x 2"+12)| + | 4] + 
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and 
pe E Za ee a) 5 ts i + (40 x 2” — 16) E + aC 
voxels] fie] 
=(0e BeBe (SESE) 


§6. Conclusion 


In this study, we have defined the first and second inverse Nirmala indices of a molecular graph. 


Furthermore, the first and second inverse Nirmala indices for certain dendrimers are computed. 
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Abstract: Let M be a noncommutative 2-torsion free semiprime I’-semiring satisfying a 
certain assumption with centre Z,(M) and T: M > M be an additive mapping. We prove 
results: 1) If T is centralizing on a Jordan [-subring J of M, then T is commutating on J; 
2) If T is centralizing right centralizer on M, then T is commutating; 3) If T is centralizing 
right centralizer on M, then T is centralizer and 4) If T is centralizing right centralizer on 
M, then T satisfies the relation 


[z, ylaBT (2) = [T(2), ylabe 


for alla,ye€ M anda,BeTl. 


Key Words: Semiprime gamma ring, semiprime gamma semiring, centralizing, left and 


right centralizers. 
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§1. Introduction 


The notion of gamma ring was first introduced in [1], which is currently notable as I'y-ring. 
Bernes [2] broadly generalized and extended the concept of I'y-ring to T-ring and shown that 
every ['y-ring is a [-ring. The [-ring is more general than the classical ring and it is concluded 
that [-ring need not to be a ring [1, 2]. Later, much theory relevant to the classical rings have 
been generalized and extended to the theory of T-rings, especially, Luh [3] and Kyuno [4] deeply 
studied on the structure of T-rings and explored various generalizations of analogous parts in 
ring theory. 

Over the years, Bell and Martindale [5] and Zalar [6] developed some notable results on 
centralizing mappings of semiprime rings. Vukman [7-10] presented may remarkable findings 
via the concept of centralizers on prime and semiprime rings. Recently, the research on cen- 
tralizers of prime and semiprime rings have been extended to prime and semiprime gamma 
rings and semiprime gamma semirings in the aspects of Jordan centralizers [12, 13], centralizers 
[11, 13C15], centralizers on Lie ideals [16, 17], centralizers with involutions [18, 19], Jordan 
derivations on Lie ideals [16] and generalized derivations on prime and semiprime gamma rings 
with centralizing and commuting [20, 21] as well. 


1Received March 3, 2021, Accepted June 8, 2021. 


42 Nondita Paul and Md Fazlul Hoque 


H.S. Vandiver introduced the algebraic study of semiring in 1934 [22, 23] and Rao [24] 
extended such research to the [-semirings and established some basic theories on gamma rings 
as well as on gamma semiring. A number of important features on I'-semirings are presented in 
(12, 25, 26]. However, the research on centralizing left/right centralizers on prime and semiprime 
gamma semiring is still unknown area. Thus the purpose of this article is to study on semiprime 
gamma semiring via centralizing right centralizers [11, 20]. The study is inspired by the work 
of [25, 26]. The results presented in this paper through out for right centralizers, which are also 
true for left centralizers because of left-right symmetry. 


§2. Preliminaries 


Let M and T be additive Abelian groups. If there exists a mapping (z,a,y) > xray of M x 
T x M > M, which satisfies the following conditions: 


(a) 
(b) ca(y + z) = xay + raz and (a + y)az = waz t+ yaz; 
(c) e(a+ B)y = way + «By; 

(d) (tay) Bz = xa(yBz) 
for all x,y,z € M anda,@ € T, then M is called a [-ring. Every ring M is a [-ring with 
M=T. 

Let M and T be two additive commutative semigroups. Then M is called a T-semiring if 
M is itself a T-ring. Obviously, every semiring M is a [-semiring with M =T[. A non-empty 
subset A of a T-semiring M is said to be a sub I-semiring of M if (A,+) is a subsemigroup 
of (M,+) and ray € A for all z,y € A anda eT. A [-semiring M is said to have a zero 
element if there exists an element 0 € M such that 0+ 4 =a =2+0 and Oar = 0 = xa for 
alla € M anda €éTI. AT-semiring M is said to be prime if ray = 0 implies x = 0 or y = 0 
for allz,y € S andaeéT. A [-semiring S is said to be semiprime if rax = 0 implies x = 0 
for alla € S andaéT. A [-semiring ™ is said to be n-torsion free if nx = 0 implies x = 0 
for alla € M. A T-semiring M is said to be commutative if cay = yax for all x,y € M and 
a€T. Let M be aT-semiring. Then the set Z7,.(M) ={xE€ M:xay=yar Vye M,aeT} 
is called the centre of the [-semiring M. Let M be aT-ring. Then [x,y], = ray — yaz is called 
the commutator of x and y with respect to a, where x,y € M andaeTl. 


We make the basic commutator identities following 


[aab, c]g = [a, c]gab + ala, B]-b + aalb, c]g, 


[a, bac|g = [a, blgac + bla, Blac + bala, e]p 
for all a,b,c € M and a, 8 € T. We consider the following assumption [11], 
aabBc = aBbac (2.1) 


for all a,b,c € M, and a,@ € T, which we extensively used in this paper. According to the 
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assumption(2.1), the above two identities reduce to 


[aab, cls = [a, clgab + aalb, cls, (2.2) 


[a, bac]g = [a, b]gac + bala, c]g. (2.3) 
The identities (2.2) and (2.3) are also used thoroughly in this article. 


An additive mapping T : M — M is called a left (right) centralizer if 
T(xay) = T(x)ay(T (way) = aT (y)) 


holds for all z,y € M anda €T. An additive mapping T : M —> M is centralizer if it is 
both a left and a right centralizer. For any fixed a € M and a €T, the mapping T(x) = aax 
is a left centralizer and T(x) = aa is a right centralizer. A mapping T : M — M is called 
centralizing if [T(x),t]a € Za(M) for alla € M, a € T. A mapping T of a T-semiring M 
into itself is said to be commuting if [T(x),z]q = 0. We recall if T : M — M is commuting, 
then [T(«), yla = [z,T(y)]a for all x,y € M and a €T. Obviously, every commuting mapping 
T :M — M is centralizing. If A be a subset of [-semiring M and ray + yax € A for all 
z,y€ Aanda€TJ, then A is called a Jordan subring of M. 


§3. Main results 


In this section, we obtain the following results with their proofs in sense of 2-torsion free 
semiprime T-semiring with the certain assumption (2.1) and using various commutation iden- 
tities. 


Theorem 3.1 Suppose M is a 2- torsion free cancellative semiprime T-semiring satisfying the 
assumption (2.1) and J is a Jordan subring of M. If an additive mapping T : M > M is 


centralizing on J, then T is commuting on J. 


Proof By the definition of centralizing T on J, we have 
(T(x), t]a € Zy(M). (3.1) 
For the linearization, we put « = x + y in (3.1), which yields 


[T(x), tla + (T(z), yla + [T(y), tla + [T(y), ya € Za(M), 
= [T(z),yla + [T(y), t]a € Za(M) (3.2) 


for allz,y € J anda€T. In particular, for y = 6x, we obtain 


(T(x), ¢Ba]y+[T(#Bxr),t]a€Z.(M) forall z,yeJ, aeTl, 
=> [T(2),t].6x + cB(T (x), cla + [T(eBx), tla € Zo(M). 
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Using the definition of the centre of I-semiring Z,(M), we have 


(T(x), a]gGa + [T (x), c]o 8a + [T(xBx),2]q € Za(M), 
=> 2(T(x),2],.6x + [T(aBx), tla € Za(M). 


Suppose « € J is a fixed element with z = [T(z),z]a € Za(M) and a 


Then (3.3) can rewrite in the following form 
[T(x),2zBa + ala =0. 
By the expansion of the commutation identities (3.4), we get 


[T'(x), 2z8a]q + [T(a), ala = 0, 

= 226[T(«), t]a + 2(T (x), 2Z]a Ba + (T(x), ala = 0, 
=> 2282+ (T(z), ala =9, 

=> (T(x), alg = —228z. 


On the other hand, we have 


[T(aBx), eB2]q € Za(M) 


for all x € J and a, €T. This implies 


[T(aBx), 2]y8a + x6[T(xBx), at]. € Za(M). 


(3.5) 


(3.6) 


(3.7) 


Now using the relation (3.7), we can write [T'(#),aBx2+xGa], = 0 and apply (3.4), it takes 


the following explicit form 


(T(x), alabx + aB8[T (x), cla + (T(x), t].6a + xB[T (x), ala = 0, 
=> —228262 + aB8z 4+ zBa+x8(—2zBz) =0, 


=> —22826x + zBat zBa—2z8z8a2 =0, by using the definition of Z,(M), 


=> 2zCa—4z28z0xa =0 > a= 228a. 


Using (3.8) in (3.5), we have 


[T(x), 2zBa]_ = —22z8z, 
=> 2{[T(z), z]aBx + zB[T(2), tla} = —22zBz, 
= (T(z), zloba + zB[T (2), tla = —2Bz, 
=> (T(z), z]oba + z6[T (2), tle = —zB2z, 

= z6[T (2), t]q = —2Bz, 


zbz = —zBz => 2zBz =0. 


(3.8) 


(3.9) 
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By the 2-torsion free semiprimeness of M, we conclude that z3z = 0 implies z = 0 for all 
8 €T. Therefore [T(x), 2]. = 0 for all  € J and hence T is commutating on J. 


Theorem 3.2 Suppose that M is a cancellative semiprime C-semiring satisfying the assumption 


(2.1). If: M — M is a centralizing right centralizer, then T is commuting. 


Proof If we consider M is 2-torsion free cancellative semiprime I-semiring satisfying the 
assumption (2.1), then the theorem is nothing to prove for J = M in account of Theorem 3.1. 
We now assume that M is not 2-torsion free -semiring. In this case, we consider the following 
relation 

2{r,T(x)lq = 0. (3.10) 


for alla € M and a €T. We now substitute x + y for x in (3.10), which yields 


2[e+y,T(t+y)la = 9, 
=> 2|2, Tyla + 2[y, T(z)le = 0, 
= [z,T(y)Jo = —ly, T(2)a (3.11) 


for all z,y € M and a €T. We now again linearise the assumption [x,T(2)]a € Za(M) by the 
transformation « = x + y, which leads to 


[z, Tyla + ly, T(®)]a € Za(M) (3.12) 
for allz,y €¢ M andae€T. By using the definition of Z,,(M) in (3.12), we enable to express as 
(Ie, Tyla + ly, T()]a), z]p = 0 (3.18) 


for all z,y € M anda, €T. This implies 


(Iz, Two + ly, T(@)Ja)Bx — 8 ([x,Ty)Ja + ly, T(2)]a) = 9, 
=> |z,Tylobe + ly, T(z)loba — c6lz,Ty)le — «Bly, T(x))q = 0. (3.14) 


Using (3.11) in (3.14), we obtain 
[z, T(y)lax + [y, T(2)]o82 + eB[x, Tyla + Bly, T(®)]a = 0. (3.15) 
Again from the assumption [z,T(x)|4 € Za(M) and the definition of Z,(M/), we found 


[x,T(t)]aBy = (waT(x) — T(x)ax) By 
= xaT (x) By — T(x)axBy 
= yBraT(x) — yBT(x)ax 
= yB(raT(x) — T(x)axr) = yblr,T(x)Ja (3.16) 
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for all z,y € M anda, ET. Applying the result(3.16) in (3.11), we get 


2[2,T (x)Ja = 0, 

=> [2,T(x)Jo + [z,T(x)Jo = 9, 

=> yBla,T(x)la + yBla, T(x) =0, right multiplying by yf, 

=> [r,T(x)|o6y + y6lx,T(x)]o =0, using (3.16), (3.17) 


for all x,y € M and a,@ €T. Adding the relations (3.15) and (3.17) and simplifying, we have 
[((wby + yb), T(x)\a + [wBx,T(y)]a = 0. (3.18) 
Now using xyy for y in the relation (3.18), we arrive at 
((eBaqy + eyyBx),T(x)Ja + [eba,T(ryy)Ja = 0 (3.19) 
for all z,y € M and a,f,7 €T. By using assumption (2.1) in (3.19), we obtain 


[(zyxBy + zyyBx),T(x)\a + [vBx,T(ryy)\a = 9, 
= [x7(aBy + yBx),T(x)la + [xbx,T(xyy)la = 9, 
=> [2,T(r)]ay(xby + yBx) + xy[(aBy + yBx),T(2)lo + 2y[xGa,T(y)Ja = 0. (3.20) 


Using (3.18) in (3.20), it reduces to 


(x)lav(xBy + yBx) = 0, 

(x)lay(xBy — yBx + 2yBx) = 0, 

(r)lay(xBy — yBx) + 2[2,T(zr)loyvyBx = 0, 

(x)Jay(aBy — yBx) =0, using (3.10), 

(x)Jarlz,y]a = 0 (3.21) 


t Yl y 
8 
HRY A NK 


for all az,y € M and a,{,7 €T. Replacing y by T(x) and 6 = a in (3.21), we obtain 
[x, T(z)la Vz, Dee) e =0 (3.22) 


for alla € M and a,y € IT. For the semiprimeness of I-semiring, [x,T(x)]ay|[z,T(x)la = 0 
implies [x,T(x)|q = 0 for all x € M and a,y € TI. Therefore T is commuting and hence the 


theorem is proved. 


Theorem 3.3 Suppose that M is a cancellative semiprime l-semiring satisfying the assumption 
(2.1). If: M — M is a centralizing right centralizer on M, then T is centralizer. 


Proof Since T is a centralizing right centralizer on M, we have 


T(xay) = caT(y) (3.23) 
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for all z,y € M anda €T. We aim to prove that T(ray) = T(x)ay for all x,y € M and 
a €T. According to the statement and Theorem 3.2, T is commuting on M. In this case, we 
can write [x,T(x)]~ = 0. This implies 


cal (x) = T(x)ax (3.24) 
for alla € M and a €T. Now linearizing the relation (3.24) by setting = x + y, yields 
[z,Tw)la + ly, T(@)]a = 0 (3.25) 


for allz,y€ M andaé€T. Applying xy for x and using the definition, we obtain 


[xBy, T(y)la + [y, T(tPy)]a = 9, 
= Bly, Tyla + [2, Ty)aby + ly, T(ey)Ja = 9, 

= [z,Ty)laSy + ly, xbT(y)la = 9, 

= xaT(y) By — T(y)axBy + yaxBT(y) — x6T(y)ay = 0, 

= xBbT(y)ay — «BT (y)ay + yaxBT(y) — T(y)axBy = 0, 

= yaxBT(y) — T(y)axBy = 0, 

= yaxBT(y) = T(y)axBy (3.26) 


for allz,y € M anda,B ET. Let T(y) = y in (3.26), then we have yaxGT(y) = yaxBy. This 
implies y8BxaT(y) = yBxay. By using the cancellation law, it shows that raT(y) = cay, which 
implies that T(xay) = xay for all x,y € M anda €T. Thus for choosing T(x) = x, we write 
T(xay) = T(x)ay. By using the assumption (2.1), this implies 


zaT(y) = T(x)ay. (3.27) 


If we consider z € M and 6 €T, then we can write (T(xay)—xaT(y))8z28(T (cay) —xaT(y)) = 
0 for all z,y € M and a € IT. By the definition of semiprime I-semiring M, it leads to 
T(xay) — caT(y) = 0. That is, 

T(xay) = xaT(y). (3.28) 


Comparing (3.27) and (3.28), we conclude that 
T(xay) = T(x)ay (3.29) 


for all z,y € M anda €TI. Therefore, T is a centralizer on M and hence, the theorem is 


proved. 


The study of the above theorems, we can provide the following remarks. 


Remark 3.1 Every centralizer on a cancellative semiprime [-semiring is commuting, because 
of T(xax) = T(x)ax = xaT (x) for all x € M and a € T and hence is a centralizing right 


centralizer. 
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Remark 3.2 An additive mapping T of a cancellative semiprime T'-semiring M is a centralizer 
if and only if it is a centralizing right centralizer on M. 


Corollary 3.1 Suppose that T is a commuting right centralizer of a semiprime T-semiring M 
satisfying the assumption (2.1), then T satisfies the relation [x, y]a8T(x) = (T(x), yla Sx for all 
z.yEeM anda,Betl. 


Proof By the statement, we have T(xay) = xaT(y) for allz,y € M anda €T. Since T is 
also commuting on M, it is easily seen that [7,T(x)], = 0. Putting « = «+ y for linearization, 
we arrive at 


[7, Tyla + ly, T(@)]a = 0 (3.30) 


for alla,y € M anda €T. Replacing y by y@z in (3.30) and using the relation [x, T(z)]q = 0, 
we obtain 


=> UYpl (ae le TT [yBa, 


) T(x)\a = 9, 
= [x, yloBT (x) + yBla, T@)a + ly, T@)abx + Bla, T(a)]a = 0, 
= [x, ylobT (x) + ly, T(x) Jax = 0, 
= (x, ylobT (x) — |T(a), yloBx = 0, 
= (2, yloBT (x) = (T(x), ylaBa (3.31) 


for all z,y € M and a,@ «IT. Hence, the theorem is proved. 


Corollary 3.2 Suppose that M is a prime T-semiring and T is a commuting right centralizer 
on M. If T(x) € Z.(M) for all x € M, then T =0 or M is commutative. 


Proof Since T(x”) € Z,,(M), in this case we have [T'(x), y]a = 0 for allz,y € M andaeTl. 
We also have [z, yla6T(x) = [T (x), yla 8x for all z,y € M and a,$ €T. Thus 


[z, y]aBT (x) =0 (3.32) 


for allaz,y € M anda,@ ET. Setting y = yyz in (3.32), we have 


[x, y¥Z]aBT(x) = 0, 

= {[z, ylavz + yz, la} BT (x) = 0, 

= (2, ylay28T (x) + yz, ZlaPT (x) = 0, 

=> [z, ylayz6T (x) =0, forall xz,y,2€ M, 
= (cay — yax)yzBT(x) = 0, 

=> (rayyz — yaxyz)bT (x) = 0. 


For the prime I’-semiringness of M, we have (tayyz — yaxyz) = 0 or BT (x) = 0. Thus we 
have seen that T = 0 or M is commutative, and hence, the theorem is proved. 
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Abstract: A robot can determine its location by a set of fixed landmarks (consider the 
robot path is like a graph). By a signal from the robot, it can determine how far it from 
each among a set of fixed landmarks. We can formulate this problem such that the robot can 
always determine its location as the following: how to compute the minimum landmarks and 
where these landmarks should be placed. The metric dimension problem of a given graph 
can answer about the above two questions. The metric basis of the graph represents the 
set of landmarks and the metric dimension of the graph is equal to the cardinality of the 
landmarks. The determining of the metric dimension of an arbitrary graph is an NP-complete 
problem. In this paper we determine the metric dimension of new classes of networks with 
polynomial algorithms. In particular, the metric dimension of the triangular snake graph, 
the ladder graph L,, Bx, the k-home chain graph Hy, and the k-kite chain graph Ky, are 
determined. We propose an exponential algorithm for finding the metric dimension of a given 
graph. We also introduce the results of computer calculations that determined the metric 
dimension of various classes of networks by using an approximate algorithm namely integer 
linear programming. Finally, a comparative analysis between the proposed approximate and 
exponential algorithms shows that the LINGO solver outperforms the proposed algorithm 
for huge graphs. 

Key Words: Metric dimension, resolving set, basis, integer programming, NP-hardness, 
robot navigation. 


AMS(2010): 05C30. 


§1. Introduction 


A robot can determine its location by a set of fixed landmarks (consider the robot path is like a 
graph). It is known that a robot is different distances from a set of fixed landmarks. By a signal 
from the robot, it can determine how far it from each among a set of fixed landmarks. We can 
formulate this problem such that the robot can always determine its location as the following: 
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how to compute the minimum landmarks and where these landmarks should be placed. The 
metric dimension problem of a given graph can answer about the above two questions. The 
metric basis of the graph represents the set of landmarks and the metric dimension of the graph 
is equal to the cardinality of the landmarks. The concept of the metric dimension has proven 
to be useful in a variety of fields. Chartrand et al. [1] applied the resolving set of a metric 
dimension in chemistry to classify chemical compounds. Khuller et al. [2] applied resolving sets 
in robotic navigation, and Sebo et al. [3] applied them in combinatorial search and optimization 
problems. 

A systematic definition of the metric dimension is proposed. Let G be a connected graph 
and d(x, y) be the distance between vertices x and y. A subset of vertices w = {wi,--: , wr} 
is called a resolving set for G if for every two distinct vertices x,y € V(G), there is a vertex 
w; € w such that d(z,w;) 4 d(y,w;). The metric dimension md(G) of G is the minimum 
cardinality of a resolving set for G. 

Figure 1 shows how the different signals determine the robot’s location. Suppose the robot 
moves on the path graph (straight line) and there is a landmark vertex (v4) that sends different 
signals (the distances are 3,2,1, and 0 between v1, v2,v3, and v4 and the landmark vertex v4, 
respectively) to the robot to determine its location. 


v4 Vo V3 V4 


Figure 1. How the different signals determine the robot’s location 


Slater [4],[5] proposed the notion of a minimum resolving set. After that Slater proposed 
a different term (location set) for the resolving set for a connected graph G. He studied the 
location number of G as the cardinality of a minimum resolving set. Harary and Melter [6] 
used a different term (metric dimension) for the same problem. The metric dimension of some 
graphs such as trees, paths, and complete graphs were determined by Chartrand et al. [7]. 
They introduced the bounds of the metric dimensions for any connected graphs. They also 
formulated this problem as an integer linear programming problem. Gerey and Johnson [8] 
proved that the metric dimension problem is an NP-complete problem for an arbitrary graph. 
The metric dimension problem for grid graphs was studied by Melter and Tomescu [9]. The 
metric dimension of graphs obtained by the cartesian product of two or more graphs was studied 
by Caceres et al. [10]. Chartrand et al. [11] find each graph of order n that has metric dimension 
ljn-—2orn-1. 

In this work we determine the metric dimension of new classes of networks with polynomial 
algorithms. In particular, the metric dimension of the triangular snake graph, the ladder graph 
Im, By , the k-home chain graph H;, and the k-kite chain graph K;, are determined. We propose 
an exponential algorithm for finding the metric dimension of a given graph. We also introduce 
the results of computer calculations that determined the metric dimension of various classes 
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of networks by using an approximate algorithm namely integer linear programming. Finally, 
a comparative analysis between the proposed approximate and exponential algorithms shows 
that the LINGO solver outperforms the proposed algorithm for huge graphs. 


§2. Main Results: Polynomial Algorithms of Some Special Classes of Networks 


Here, we prove that the metric dimension of the triangular snake graph A; — snake, the ladder 
graph L,, the k-home chain graph Hx, the k-kite chain graph K, and the k-envelop chain 
graph E,, equals 2. Samir Khuller et al [12] proved that for graph G with metric dimension 2 
and metric basis {a, b}, the following are true: 

(1) There is no more than one path (P) between a and }; 

(2) The degree (a) < 3 and degree (b) < 3; 

(3) The degree (c) < 5 such that ce P,c#AaandcH b. 


Theorem 2.1 Let G_ be a triangular snake graph (A, — snake) with k blocks and n vertices, 
then md(G) = 2. 


Figure 2. Triangular snake graph (A; — snake). 


Proof We label the triangular snake network G = (A, — snake) as shown in Figure 2 such 
that & is the blocks number. It is clear that | V(G) | isn = 2k +1. Let w = {v1, vox} so that 
the proof has three cases, i.e., k = 1, k > 1 with the odd labeling of the vertices and k > 1 with 
the even labeling of the vertices. 


Case 1. For k = 1, the proof is trivial because the graph G = C3. 


Case 2. For k > 1, the odd labeling of the vertices is the following: 
Begin 
for (i =1,i<n—2;i=i+2) do 
d(v;,w) = (5+, 3) 


2 2 


end 


(vpn, w) = (S21) 
End 
Case 3. For k > 1, the even labeling of the vertices is the following: 
Begin 
for(i = 2;i <n—3;i=i+2) do 
a(vy,w) = (§, $4) 
end 


d(vpn_-1,W) = (S52,0) 
End 
This completes the proof. 
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Obviously, there are no two vertices with the same labeling , we then obtain a resolving 
set w with | w|, so we have md (A, — snake) = 2. The proof of the algorithm of Theorem 
1 involves a for-loop, so the algorithm complexity is O(n), indicating that it is of polynomial 


time. 


Theorem 2.2 Let G be a ladder graph L,, where n > 4, then md(G) = 2. 
Y, Y, 
y, ¥. V. bale 2 
1 2 3 — 1 - 
Va Va Veo "ais ms 
2 2 
Figure 3. Ladder graph L,,. 


Proof We label the triangular snake network G = (A; — snake) as shown in Figure 2 such 
that & is the blocks number. It is clear that | V(G) | isn = 2k +1. Let w = {v1, vox} so that 
the proof has three cases, i.e., k = 1, k > 1 with the odd labeling of the vertices and k > 1 with 
the even labeling of the vertices. 


Case 1. For k = 1, the proof is trivial because the graph G = C3. 


Case 2. For k > 1, the odd labeling of the vertices is the following: 
Begin 
for(i = 1;i <n-—2;i=i+2) do 


d(vi,w) = (#5) 


d(vn,w) = (22,1) 


Case 3. For k > 1, the even labeling of the vertices is the following: 
Begin 
for(i = 2;i <n—3;i=i+2) do 
d(vy,) = (, $2) 
end 


d(Un—1,w) = (“52,0) 
End 
This completes the proof. 


We label a ladder graph G = L,, as shown in Figure 3. It is clear that | V(G) | isn = 2k+2 
such that k is the blocks number of G. Let w = {v1, vg+i} 
Begin 
for (@=1i< 


end 
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for (i =n;i < $+1;i—-) do 
d(u;,w) = (w—4+1,4- 8) 
end 
End 
Obviously, there are no two vertices with the same labelling, we then obtain a resolving 
set w with | w|, so we have md(L,,) = 2. The algorithm of the proof of Theorem 2 contains two 
for-loops, but they are not inner loops, so the algorithm complexity is O(n), indicating that it 
is of polynomial time. 


Theorem 2.3 Let G be graph By, , where n > 4, then md(G) = 2. 


1 Vv 
nn. ” 
4“ %s i oz 
/|T{[ itil 
e 
V,, V,, a] 


Figure 4. Graph B, 


Proof We label the triangular snake network G = (A; — snake) as shown in Figure 2 such 
that & is the blocks number. It is clear that | V(G) | isn = 2k +1. Let w = {v1, vox} so that 
the proof has three cases, i.e., k = 1, k > 1 with the odd labeling of the vertices and k > 1 with 
the even labeling of the vertices. 


Case 1. For k = 1, the proof is trivial because the graph G = C3. 


Case 2. For k > 1, the odd labeling of the vertices is the following: 
Begin 
for (i =1,i<n—2;i=i+2) do 
d(vi,w) = (, *5*) 


end 


(vy, w) = (22.1) 
End 


Case 3. For k > 1, the even labeling of the vertices is the following: 
Begin 
for (i = 2;i<n—3;i1=i+2) do 
d(vi,w) = (5,3) 
end = d(vn_1,w) = (252.0) 
End 
This completes the proof. 


We label a graph G = B, as shown in Figure 4 such that k is the blocks number. It is 
clear that V(G) isn = 2k +4. Let w = {v1, vg42}. We want to prove that w is a resolving set 
by showing that d(v,,w) 4 d(v;, w) for all i 4 7. Observe that: 

Begin j; =0; j2=1 
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for (¢=1,i<k+1;i++) do 
A(v;,w) = (Hrs J2) 
g0 =i Re = Je oad 
end 
Uver2, w)= (1,0) jr =15 j2 =1 
for (@=k+1;i<n—1;i++4+) do d(v;,w) = (j1,J2) 
ji = Jr rlyje= jo +l 
end 
d(vn,w) = (737, 3) 
End 
Obviously, there are no two vertices with the same labelling, we then obtain a resolving 
set w with] w|, so we have md(B,;) = 2. The algorithm of the proof of Theorem 3 has two 
for-loops, but they are not inner loops, so the algorithm complexity is O(n), indicating that it 
is of polynomial time. 


Theorem 2.4 Let G be a k-home chain graph Hj, , where k > 2 then md(G) = 2. 


Vz V2 Vg 


V2ke4 


Voke2 Seer 
Vanes Vaxea 311 ake 


Figure 5. k-home chain graph Hx. 


Proof We label the triangular snake network G = (A; — snake) as shown in Figure 2 such 
that & is the blocks number. It is clear that | V(G) | isn = 2k +1. Let w = {v1, vox} so that 
the proof has three cases, i.e., k = 1, k > 1 with the odd labeling of the vertices and k > 1 with 
the even labeling of the vertices. 


Case 1. For k = 1, the proof is trivial because the graph G = C3. 


Case 2. For k > 1, the odd labeling of the vertices is the following: 
Begin 
for (¢ = 1l,i<n—2;i=i+4+2) do 


Case 3: for k > 1, the even labeling of the vertices is the following: 
Begin 
for (i = 2;i<n—3;i=i+2) do 
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d(u;,w) = (3, =") 
end = d(vn_1,w) = (252.0) 
End 
This completes the proof. 


We label a graph G = Hy, as shown in Figure 5 such that k is the blocks number. It is 
clear that | V(G) | isn = 3k + 2. Let w = {v1, ux} .We want to prove that w is a resolving set 
by showing that d(v,,w) # d(v;, w) for all i 4 7. Observe that: 

Begin d(v,,W) = (0,k) j=k-1 

for (¢ = 2;i<k—1;i++) do 
d(vi,W) = (4,9),3 =F-1 
end 
d(ve, W) = (k,0), d(veszi, W) = (1, k) 
for (i =k+2;i < 2k;i++) do 
fi=l1jje=k-1 
d(vi, W) = Ga, J2) 
Jiv=ji+13j2=J2 -1 
end 
d(vop41,W) = (k,1) d(von42,W) = (2,k +1) 
for (¢ = 2k+3;i<n—1;1++) do 
ji =2,j2=k 
d(vi, W) = Gi, Je) 
ji=jitl1jj2=Jj2 -1 
end 
d(vm,W) = (£4, 2) 

End 

Obviously, there are no two vertices with the same labelling; we then obtain a resolving 
set w with | w|, so we have md(H;,) = 2. The algorithm of the proof of Theorem 4 has three 
for-loops, but they are not inner loops, so the algorithm complexity is O(n), indicating that it 
is of polynomial time algorithm. 


Theorem 2.5 Let G be a k-kite chain graph Ky, where k > 2 with n vertices and k blocks, then 
md(G) = 2. 


p 


, ) , 
2n+2 J 2n+5 1 2n+8 Vis n 


3 3 3 


Figure 6. k-kite chain graph K,. 
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Proof We label the triangular snake network G = (A; — snake) as shown in Figure 2 such 
that & is the blocks number. It is clear that | V(G) | isn = 2k +1. Let w = {v1, vox} so that 
the proof has three cases, i.e., k = 1, k > 1 with the odd labeling of the vertices and k > 1 with 
the even labeling of the vertices. 


Case 1. For k = 1, the proof is trivial because the graph G = C3. 


Case 2. For k > 1, the odd labeling of the vertices is the following: 
Begin 
for (¢ = 1,i<n—2;i=i+2) do 


Case 3: for k > 1, the even labeling of the vertices is the following: 
Begin 
for (i = 2;i<n—3;i1=i+2) do 
atest) = (§25$4) 
end 
A(Un-1,W) = (“52,0) End 


This completes the proof. 


We label the k-kite chain graph G = Ky, as shown in Figure 6 such that k is the blocks 
number of Ky. It is clear that | V(G) | isn = 3k +2. Let w = {v1, vp4i} . We want to prove 
that w is a resolving set by showing that d(v;,w) 4 d(v;,w) for all #7. Observe that: 

Begin 

for (¢=1;i<k+1;i++) do 

d(u,,W) = (¢-1,k -—i+1) 
end 

fi=lje=k 

for (t=k+2;i <2k+1;i++) do 

d(u;,W) = (f1, 92) 31 =J1 +15 2 =Jj2-1 
end 

jx =ljja=k+1 
for (¢ = 2k + 2;1 < 3k+1;i++) do 
d(vi,W) = (93,54); Js =Js t+ lija=js2—-1 

end 

End 

Obviously, there are no two vertices with the same labelling, and so we obtain a resolving 
set w with| w|. Thus,md(k;,) = 2. The algorithm of the proof of Theorem 5 has three for- 
loops, but they are not inner loops, so the algorithm complexity is O(n), indicating that it is 
of polynomial time. 
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§3. An Exponential Algorithm for a Given Network 


In the previous section, we proposed polynomial time algorithms for special cases (triangular 
snake graph A; — snake, the ladder graph L,, By, the k-home chain graph Hy, and the k- 
kite chain graph K;,). Here, we introduce an algorithm that gives the metric dimension of an 
arbitrary graph G. Unfortunately, the time complexity of this algorithm is exponential time 
O(2”). Recall that Gerey and Johnson [8] showed that determining the metric dimension of an 
arbitrary graph is an NP-complete problem. 


Algorithm 1. Finding the metric dimension of a given graph 


Input : An adjacency matrix A[n][n] of an n-vertex simple connected graph G. 
Output: The metric dimension of G. 
Begin Apply Floyed -Warshall’s method to compute the distance matrix D[n][n] of G. 
Initialization: 
Sien = 0; Enan = 0; counter = cardinality = 0 and metric_dimesion = inf 
Generating all subsets 
while ( counter < 2” — 1) 
for (i= n;i > 1;i—) do 
if S[i] =0 then 
Sit] =1 
maz =i+tl 
for (j = max;j <n;j7 ++) do 
S[j] = 0 
end 
counter = counter +1 
cardinality = non-zero elements of S[i| 
break; 
end 
end 
Check all subsets of P(Vgq)whether resolving set or not. 
for (i =n;i > 1;i—) do 
for (j =; j 2 1;j-) do 
Efil[, 7] = 9 
end 
end 
for (i =n;i > 1;i—) do 
for (j =n;j = 1;j—) do 
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if S{j] = 1 then 


end 
for (i =n;i > 1;i—) do 
for (j =i+1;j >1;j—) do 
if E/i][:] = Elyl[:] then 
break 


end 
end 
if ( E[é][:] or Ely][:] ) and (metric_dimension > cardinality) then 
metric_dimension = cardinality; 
end 
end while 


End Begin 


Example 3.1 We show the intermediate stages of above Algorithm for the cycle graph. It is 
clear that the graph G = C’3 has 3 vertices and 3 edges. The graph G has adjacency matrix 


011 
A and distance matrix D such that A= D=]1 0 1 |. Algorithm 1 can calculate 2” 
1 1 0 


subsets as follows: 
S=[0 0 0]%=[0 0 1]%=[01 0]%=[01 1] 


Si=([1 0 0]S=[1 0 1]Se=[1 1 o]S=[1 1 1] 


Now, checks each of the above subsets S$; : 7 = 0,--- ,2” — 1 to determine whether a 


resolving set as follows: the matrix E; is calculated for each subset $; such that: 
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If there exist two rows are equal in the matrix E£; then the subset $; is non-resolving set 
otherwise it is resolving set; the algorithm determines that $3, 55,5‘ and S7 are resolving sets 
and the others are non-resolving sets. Finally Algorithm 1 determines the metric dimension 
that is equal to the number of elements of the smallest resolving set for G so md(G) =| $3 |=| 
55 |=| Se |= 2. 


Complexity of Algorithm 1. Obviously, Algorithm 1 is an exponential algorithm. It consists 
of one while-loop that has four for-loops each with an inner loop. So the total complexity of 
Algorithm 1 is 


§4. Formulation of the Problem as Integer Linear Programming Model 


In the previous section, we proposed an exponential algorithm that determines the metric 
dimension for a given graph but this algorithm cannot determine the metric dimension of very 
large graphs in a reasonable time. In this section we introduce a powerful technique “integer 
linear programming technique” that determines the metric dimension for a very large graph. 
This technique finds the metric dimension in a reasonable time. We now describe this problem 
of determining the metric dimension and a basis for a graph in terms of an integer linear 


programming problem [13-16]. Chartrand et al [11] formulated the problem of finding the metric 


dimension as follows: Let G be a connected graph of order n with V(G) = {v1, v2,--: , Un} and 
let D = [d;;] be the distance matrix of G, that is, dj = d(vi,v;) for 1 < i337 <n. For 
x, € {0,1}; 1 < i < n, define the function z by: z(a1,29,--- ,%n) = >>} 2; such that the 


number of constraints equals to ;then, the integer linear programming will be as follows: 


n! 
2!(n—2)! 
Minimum z= 27, +%24+-:-+2n 


Subject to: 
for (i = 1,i=n; i++) do 
for (j = 1, j=n; j++) do 
| diy — dj | 21+ | dig — djg | vo+ | diz — dj3 | 3 > 0 
end 


end 


where x; € {0,1}. 
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Example 4.1 From the above formulation, we have the following mathematical model for 
the cycle graph C3 with three vertices. min z = 2, + %2+ 2x3 subject to 02; + 2g + 23 > 0, 
21+ 0a +23 > 0, 1 + 22 +023 > 0 and 21,22,23 € {0,1}. When we solve the above 
mathematical model, we obtain the objective function z = 2 so the metric dimension for C3 is 
2. 


85. Numerical Experiments 


A connected graph which consists of k- blocks (triangles) and the block-cut-point graph is a 
path is called k-triangular snake (or A,-snake) [12]. A connected graph in which the k& blocks 
are isomorlphic to the cycle C’, and the block-cutpoint graph is a path is denoted by kC,,- 
snake [13]. The ladder graph L,, is the Cartesian product of P; and P, is the path graph 
on n nodes. We describe our numerical experiments and present the computational results, 
which demonstrate the efficiency of the proposed algorithm on a set of test problems (path, 
cycle, ladder and A;-snake). Table 1 describes the computing environment. MATLAB solver 
was used to solve the mathematical model. 


Table 1. Description of the computing environment 
CPU Intel (R) Core (TM) i3-3217U CPU@ 1.80 GHz 
RAM Size 4 GB RAM 
MATLAB version R2018a (9.4.0.813654) 


Table 2 and Figure 6 show that the LINGO solver outperforms the proposed exponential 
Algorithm 1, with respect to CPU time, for determining the metric dimension for the cycle 
and path graphs. The LINGO solver can solve the problem for graphs with large sizes in a 
reasonable time. The proposed exponential Algorithm 1 can find the metric dimension for a 
given graph but in a non-reasonable time. However, the LINGO solver is not guaranteed to 
obtain the exact metric dimension for general arbitrary graph so it is an approximate algorithm. 


LINGO solver against proposed algorithm for cycle graph LINGO solver against proposed algorithm for path graph 
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Figure 7. A comparison between the proposed algorithm and LINGO solver for finding the 
metric dimension for cycle and path graphs 
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Table 2. A comparison between the proposed algorithm and LINGO solver for finding the 


metric dimension for cycle and path graphs 


Cycle Path 
d | md | uneo cpu | Aigorithmi cpu | d | md | uinco cpu | Agorithm cpu 
1} —-] - - 0 1 0.00102 0.00173 
2);-|- - 1 1 0.00103 0.00436 
3 | 1 2 0.03165 0.0060 2 1 0.02167 0.00573 
4 | 2 2 0.03816 0.00602 3 1 0.02818 0.00584 
5 | 2 2 0.04210 0.00621 4 1 0.03212 0.00781 
6 | 3 2 0.04324 0.00817 5 1 0.03326 0.00861 
7 | 3 2 0.04405 0.01283 6 1 0.03408 0.01692 
8 | 4 2 0.04956 0.02309 7 1 0.03959 0.02898 
9 | 4 2 0.04966 0.04872 8 1 0.04844 0.05044 
10 | 5 2 0.04977 0.11682 9 1 0.04975 0.12025 
11] 5 2 0.04979 0.27877 10) 1 0.04979 0.27940 
12 | 6 2 0.04985 0.76805 11) 1 0.04985 0.66518 
13 | 6 2 0.04999 2.01020 12) 1 0.04996 1.55260 
14) 7 2 0.05397 3.69079 13) 1 0.05386 3.64734 
15 | 7 2 0.05429 8.48939 14) 1 0.05438 8.45053 
16 | 8 2 0.05756 18.02579 15) 1 0.05764 20.24835 
17/ 8 2 0.05981 44.24297 16) 1 0.05993 45.47590 
18 | 9 2 0.06087 108.54759 | 17] 1 0.06096 102.98411 
19 | 9 2 0.06445 241.59339 | 18) 1 0.06455 251.55695 
20 | 10} 2 0.06762 511.26855 | 19} 1 0.06776 569.33483 
soon LINGO solver against proposed algorithm for -snakes Graph 
oil sap , LINGO solver against proposed algorithm for ladder graph 
200 f Proposed Ag 


Number of blockes 


Number of blockes 


10 


Figure 8. A comparison between the proposed algorithm and LINGO solver for finding the 


metric dimension for A;,-snakes and ladder graphs 
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It is clear that the curve for the LINGO solver is not visible in Figure 7, Figure 8 and 
Figure 9 because it runs essentially along on the z-axis. Table 3 and Figure 8 show that the 
LINGO solver outperforms the proposed algorithm, with respect to CPU time, for finding the 
metric dimension for the A;-snake and ladder graphs. 


Table 3. A comparison between the proposed algorithm and LINGO solver for finding the 
metric dimension for A;-snakes and ladder graphs 


Ag-snakes graph Ladder graph 

k | n | d | md | unco cpu | aigorithmi cpu | k | nm | d | md |} unco cpu | algorithm cpu 
1} 341 2 0.05337 0.00674 1 | 4 | 2 2 00.2791 0.00930 
2/5 | 2 2 0.05624 0.00718 2 | 6) 3 2 0.03037 0.01330 

3°] 7 | 3 2 0.05803 0.01141 3.| 8 | 4 2 0.03301 0.02808 
4/9 ]4 2 0.05951 0.03922 4 | 10) 5 2 0.03336 0.11058 

5 | 11] 5 2 0.06009 0.20711 5 | 12) 6 2 0.03389 0.62880 

6 | 13] 6 2 0.06476 1.19454 6 | 14) 7 2 0.03431 3.71727 

7 | 15) 7 2 0.06445 6.00546 7 | 16] 8 2 0.03594 19.96398 
8 | 17] 8 2 0.06501 31.78448 8 | 18) 9 2 0.04105 107.23148 
9 | 19] 9 2 0.06679 167.52252 9 | 20); 10} 2 0.04293 651.14236 
10 | 21} 10] 2 0.07332 994.01483 | 10 | 22] 11) 2 0.05581 | 6233.951601 


Table 4 and Figure 9 show that the comparison between the proposed exponential Algo- 
rithm 1 and the LINGO solver for determining the metric dimension for complete graphs with 
a metric dimension larger than 2. The LINGO solver outperforms the proposed exponential 
Algorithm 1 in terms of the CPU time. There is no guarantee of the superiority of the LINGO 
solver in determining the exact metric dimension for a general arbitrary graph. The proposed 


exponential Algorithm 1 has this guarantee for a given graph. 


LINGO solver against the proposed algorithm for complete graph 
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Figure 9: A comparison between the proposed algorithm and LINGO solver for finding the 
metric dimension for A,-snakes and ladder graphs 
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Table 4. Comparison between the proposed algorithm and LINGO solver for finding the 
metric dimension for complete graphs. 


complete 
n | d | md | LINGO CPU | Algorithm1 CPU 
1 WO) 4 0.00101 0.001000 
Be | inst 0.00102 0.004620 
aml ae 0.02165 0.004850 
Ac \ ids | 3 0.02815 0.006700 
5/1] 4 0.03210 0.009000 
Gelli Ni a5 0.03320 0.016920 
7/1] 6 0.03401 0.029990 
Sag 0.03952 0.050440 
Gel ds 0.04820 0.120250 
16>] ||| “9 0.04943 0.279400 
1h | d'l\| 10 0.04968 0.665180 
9) 1) ||| ot 0.04992 1.552600 
13.) 4 || 48 0.04998 3.140139 
TA. |Pt-\|| 43 0.05385 3.647340 
15 | 1] 14 0.05447 3.819237 
16/1) 15 0.05782 5.596174 
17/1] 16 0.05905 10.697527 
18) 1] 17 0.06079 15.722732 
19 | 1] 18 0.06444 17.793874 
20| 1] 19 0.06782 67.070074 


§6. Conclusion 


We determined the metric dimension of new classes of networks with polynomial algorithms. 
In particular, we proved that the metric dimension of the triangular snake graph. A; — snake, 
the ladder graph L,, By, the k-home chain graph Hy, and the k-kite chain graph Ky, equals 
2. We proposed an exponential algorithm for finding the metric dimension of a given graph. 
We also introduced the results of computer calculations that determined the metric dimension 
of various classes of networks by using an approximate algorithm, namely, integer linear pro- 
gramming. Finally, a comparative analysis between the approximate algorithm (the LINGO 
solver) and the proposed exponential algorithm showed that the LINGO solver outperforms the 
proposed algorithm for very large graphs. In future work, we will apply variant integer linear 
programming models [23-24] on the variant networks [25-26]. 
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Abstract: A graph is k-degenerate if its vertices can be successively deleted so that when 
deleted, they have degree at most k. A k-tree is a graph that can be formed by starting with 
Kr+1 and iterating the operation of making a new vertex adjacent to all the vertices of a 
k-clique of the existing graph. A structural characterization of maximal 2-degenerate graphs 
with diameter 2, containing 45 distinct infinite classes of graphs, is proven. A forbidden 


subgraph characterization of k-trees with diameter 2 is proven. 
Key Words: Degeneracy, diameter, k-tree, k-path. 
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§1. Introduction 


In this paper, we work toward a characterization of the maximal k-degenerate graphs with 
diameter 2. 


Definition 1.1 Let k be a positive integer. A graph is k-degenerate if its vertices can be 
successively deleted so that when deleted, they have degree at most k. A graph is maximal 
k-degenerate if no edges can be added without violating this condition. 

A k-tree is a graph that can be formed by starting with Ky41 and iterating the operation of 
making a new vertex adjacent to all the vertices of a k-clique of the existing graph. 


A k-leaf is a degree k vertex of a maximal k-degenerate graph. 


Lick and White introduced k-degenerate graphs in 1970 [13], and their properties have 
been studied by many authors [2, 7, 8, 9, 10, 11, 12, 14, 16, 19]. For n > k+1, a maximal 
k-degenerate graph has at least one k-leaf, and a k-tree has at least 2. 

The three maximal 2-degenerate graphs of order 5 are shown below [3]. The two on the 
left are 2-trees. 


Kn + Ks Pyt+ ky Ws 
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Undefined notation and terminology will generally follow [3]. In particular, the join of 
graphs G and H is denoted G+ H, and the distance between vertices u and v is d(u,v). The 
eccentricity eg (v) of a vertex v is the maximum distance between v and any other vertex of 
G. If G is a graph, the square G? is formed by adding all edges between pairs of vertices with 
distance 2 in G. 

We solve two special cases of the problem of characterizing the maximal k-degenerate 
graphs with diameter 2. One restricts the problem to maximal 2-degenerate graphs, the other 
restricts it to k-trees (which are all maximal k-degenerate). The first provides a structural 
characterization, and the latter provides a forbidden subgraph characterization. 

This work is inspired by a previous paper [6]. I coauthored with Zhongyuan Che on the 
Wiener index of maximal k-degenerate graphs. We showed that the Wiener index is minimized 
when these graphs have diameter 2. We also characterized 2-trees with diameter at most 2. 


Proposition 1.2([6]) The following are equivalent for a 2-tree G: 


(1) G has diameter at most 2; 

(2) G does not contain P?; 

(3) G is T+ Ky for any tree T, or any graph formed by adding any number of vertices 
adjacent to pairs of vertices of K3. 


§2. Maximal 2-Degenerate Graphs with Diameter 2 


In this section, we provide a structural characterization of maximal 2-degenerate graphs with 
diameter 2. 


Definition 2.1 A dominating vertex of a graph is a vertex adjacent to all other vertices. A 
fan is the graph P,-1 + Ky. 


Lemma 2.2 If G is a maximal 2-degenerate graph with order n > 3 containing a dominating 
verter, then G is a 2-tree that can be represented as T + Kk, for some tree T. If G has exactly 


two 2-leaves, then it is a fan. 


Proof We use induction on n. When n = 3, G = Kz and the result holds. Let G be 
a maximal 2-degenerate graph with order n containing dominating vertex u, and assume the 
result holds for all graphs with order n—1. Then G has a 2-leaf v, which is adjacent to u. Now 
G —v is maximal 2-degenerate with order n — 1 [13], so it is a 2-tree that can be represented as 
T+ Kk ,. Then the other neighbor of v is a neighbor of u, so G is a 2-tree that can be represented 
as T+ Ky. 

If G has exactly two 2-leaves, then deleting its dominating vertex produces a tree with 


exactly two leaves, a path. Thus G is a fan. 


Definition 2.3. When constructing a maximal 2-degenerate graph, we duplicate a 2-leaf by 
adding another 2-leaf with the same neighborhood. The inside graph H of a mazimal 2- 
degenerate graph G is formed by deleting all the 2-leaves. The stem set of G is the set of 
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neighbors of 2-leaves. 


Note that in a maximal 2-degenerate graph with diameter 2, any 2-leaf can be duplicated 
arbitrarily many times. The new 2-leaf is distance two from its duplicate, and hence at most 
two from every other vertex. Thus the result is a maximal 2-degenerate graph with diameter 2. 


Lemma 2.4 In any maximal 2-degenerate graph with diameter 2 and order n > 3, either 


(A) all 2-leaves have a single common neighbor, or 
B) the stem set is S = {u,v,w}, and there are 2-leaves with neighborhoods {u,v}, {u, w}, 
and {v, wh}. 


Proof Any maximal 2-degenerate graph with diameter 2 has at least one 2-leaf. No 2- 
leaves can have disjoint neighborhoods, since then they would be at least distance 3 apart. 
If all 2-leaves have the same neighborhood, the result follows. If two 2-leaves have distinct 
neighborhoods, we may call them {a,b} and {a,c}. Any other 2-leaf must have neighborhood 


{b,c} or {a,x} for some zx. 


Theorem 2.5 Let G be a maximal 2-degenerate graph with diameter 2. Then G is a 2-tree 
that can be represented as T + Kk, for some tree T, or the inside graph of G is one of the 44 
possibilities shown below. (Vertices labeled x may be duplicated arbitrarily many times.) There 
must be at least one 2-leaf of G neighboring any pair of black vertices or pair of black and gray 
vertices, and there may be at least one 2-leaf of G neighboring any pair of black and lightgray 


><> <— 
<N <N <u <x 
<1 <i << <2 
ro 


vertices. 


— 
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The proof of this theorem has many cases. We use Case A.2.1 to mean case A, Subcase 
2, Subsubcase 1, and similarly for the other cases. Figures are referenced in parentheses, with 
labels beginning with their main case (A or B). We say an inside graph is valid if it is the inside 
graph of a maximal 2-degenerate graph with diameter 2. 


Proof Let G be a maximal 2-degenerate graph with diameter 2 with inside graph H. By 
Lemma 2.4, there are two possibilities for the positions of the 2-leaves. 


Case A. All 2-leaves of G have a single common neighbor uw. 


Case A.1 If u is a dominating vertex of H, it does the same for G, so by Lemma 2.2, G 
is a 2-tree that can be represented as T+ kK, for some tree T. 


Case A.2 If wu has eccentricity 2 in H, let v1, ... v; be distance 1 from u, wy, --- , we 
be distance 2 from u. Now no 2-leaf of H has neighborhood {u,v;} since a 2-leaf of G that 


neighbors it and u is more than 2 from w}. 


Case A.2.1 If w, is a 2-leaf of H, there is a 2-leaf of G that neighbors it and u. Then w 
neighbors all other w;, and since w; neighbors some v;, k < 2. If k = 1, then u is a dominating 


vertex of H — w,. By Lemma 2.2, H — wy is a 2-tree. Now its 2-leaves are not 2-leaves of H, 
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aside from possibly u. Then wy is adjacent to all (two) of them, and H — wy is a fan with at 
most five vertices (Al, A2, A3). 

Since all 2-leaves of G have a single common neighbor u, it is colored black (uniquely, in 
Case A). Any 2-leaf of H must be black or gray, and any vertex distance 3 from wu will be gray. 
If {u, u’} is a dominating set of H, then u’ will be lightgray if not already colored. Since these 
statuses are trivial to check, verification will be left to the reader for the other figures. 


Al A2 A3 


If k = 2, there is no 2-leaf of H with neighborhood {u, w2}, since a 2-leaf of G neighboring it is 
not within 2 of w;. Then we is a 2-leaf of H — w,. As before, H — w, — we is a 2-tree, and w 1 
and wz have two or three neighbors in it, including all its 2-leaves. Now T = H — w, — wg —u 
is a tree with all vertices either neighbors of w2 or within 2 of wy. 

If T a path, its length is at most 5. If T = Pp, there is one possibility (A4). If T = P3, 
w, may neighbor a leaf and w2 may or may not neighbor the nonleaf, or w,; may neighbor the 
nonleaf (A5, A6, A7). If T = Py, w; may neighbor a leaf or nonleaf (A8, A9). If T = Ps, wi 
must neighbor the middle vertex, and we neighbors the leaves (A10). If T has three leaves, w2 
neighbors two, and w) neighbors the third, so T = Ky,3 (A11). 


NN << 
ae a 


Case A.2.2 Suppose there is a 2-leaf v; of H neighboring u and w;. Then there is a 
2-leaf of G neighboring u and v;. Then there is no wa, but v; may be duplicated arbitrarily 
many times. Let K be the inside graph of H (delete v, and all its duplicates). Then wy, is a 
2-leaf of K. Then wu is a dominating vertex of kK — w 1, so by Lemma 2.2, K — wy is a fan. This 
fan must have order 3, 4, or 5 (A12, A13, Al4). 


po <> 


Case A.2.3 If u is a 2-leaf of H and no w; is, 7 = 2. If both v; and v2 are 2-leaves of 
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H —u, then H — u— v, — v2, has order at most 4, so it is Ky (A15), K3 (A16), or K4—e. In 
the latter case, there are two ways to attach v, and v2 to K,—e (A17, A18). 


A ot <A Hw 


Assume v, is a 2-leaf of H — u and vg is not. If vy © v2, say wy © v1. Then v2 is adjacent to 
all other w’s. If vo © wy}, v2 is adjacent to all vertices, so by Lemma 2.2, H is a 2-tree, and 
some w; is a 2-leaf, contrary to assumption. If vg wy}, then wy, is a 2-leaf of H —u—v,. By 
Lemma 2.2, H —u— v1, — wy, is a fan. Now some 2-leaf of G has neighborhood {u, w;}, so all 
ws must be adjacent, and & = 3 (A19). 

Assume v1 v2. Since v1 is a 2-leaf of H — u, its neighbors are (say) w1 and w2. Now v2 
is adjacent to all other w’s, and k > 2. Now some 2-leaf of G has either v2 or w; as a neighbor, 
so one of these vertices neighbors all w’s (excluding itself). Then H — u—v has a dominating 
vertex, so by Lemma 2.2, it is a fan with 2-leaves w, and we. If v2 is the dominating vertex, 
the fan has order at most 5, due to v;. Order 5 duplicates A14, but order 4 yields a new case 
(A20). If (say) w3 is the dominating vertex, the fan has order 5 or 6 (A21, A22). 


<x << <X 


Case A.3 If ey (u) > 2 and vertex y is at least 3 from u, then {u, y} is the neighborhood 
of a 2-leaf a of G. If dy (u,y) > 4, there is a vertex z with dy (u, z) = 2 and dy (a, z) > 2, so 
this is impossible. Thus ey (u) = 3. Let v1, ... vj; be distance 1 from u, wi, ... we be distance 
2 from u, and 7, ... 2, be distance 3 from u. Note j,k > 2 since H has no cut-vertex [13]. 

Now all vertices in the stem set other than u must be adjacent to each w; and 2; (else a 
2-leaf has eccentricity more than 2). No v; is in the stem set, since it cannot be adjacent to an 
x;. Since Ky, is not 2-degenerate, there are at most 3 stems excluding u, and ! < 2. No w; is 
a 2-leaf of H, since if there were, it would be adjacent to a v;, and all w; and z;. Now 2, isa 
2-leaf only if there is no x2, so H has at most two 2-leaves. 


Case A.3.1 Assume wu and 2x; are 2-leaves of H. Then j = k = 2, and there is no 23. 
Thus H has order 6, and H —u— 2, = K4—e. There are three ways it can be arranged, but 


the case where w, «* wz combines into the case where v; vg. In the third case, H = Pe 


(A23, A24). 
A23 A24 
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Case A.3.2 Assume u is the only 2-leaf of H, and | = 1 (there is no x2). Then at least 
one of v; and v2 are 2-leaves of H — u. If both are 2-leaves, then 3 < k < 4 since each w, is 
adjacent to some v;. If k = 3, then H-—u—v,—v2 = K4—e by Lemma 2.2. Then v; and v2 have 
one common neighbor, and there are two choices (A25, A26). If k = 4, then H — u— v1 — vo is 
Py+ ky or ky 3+ 4h, by Lemma 2.2. If it is Py+ 44, there are three choices for the adjacencies 
between the v’s and w’s, two of which produce valid inside graphs (A27, A28). If it is Ky3+F4, 
some v and w have distance more than 2. 


A25 A26 A27 A28 


Assume only v; is a 2-leaf of H — u. If its neighbors are v2 and (say) wi, at least one of which 

are 2-leaves of H — u— v1. If vg is a 2-leaf of H — u— v1, k = 3, and its neighbors are either 

adjacent or not (A29, A30). If ve is a not 2-leaf of H — u— v1, wy, is, with neighbors x; and 

vg or (say) we. If wy ve, 1 and ve are adjacent to all remaining w’s. Thus wg is the only 

2-leaf of this graph, which is W; (A31). If wi we, x1 and v2 are adjacent to all w’s of 
—u—v,— wy. Thus wy is the only 2-leaf of this graph, which is W; (A32). 


(SSNS SS 


Suppose v is the only 2-leaf of H — u with neighbors (say) w; and we, and wy is a 2-leaf of 
AH —u-—vy,. If w; has neighbors x; and v2, then H — u— v, — wy has order at least 4. Now w2 
is adjacent to all other w’s (so v; is distance 2 from them) and v2 is adjacent to all w’s, except 
perhaps we. Since x, is adjacent to all w’s, H — u—v, — w, contains K3,2, so k < 4. There 
are two possibilities (A33, A34). If w; has neighbors w3 and 21, then ws neighbors v2 and 2}. 
As before, H — u — v1 — wi — wz contains K3~,~3, so k <5. There are two possibilities (A35, 


A36). 
A33 A34 A35 x 


Case A.3.3 Assume u is the only 2-leaf of H, and 1 = 2. Then 2 < k < 4. Now one 
or both of vy and vz are 2-leaves of H — u. If k = 2, there are two cases, both leading to 
valid graphs (A37, A38). If k = 3, there is one way to make both v, and v2 2-leaves of H — u. 
However, some v and w will have distance more than 2, so this is not to a valid graph. If only 
v1 is a 2-leaf this leads to a valid graph (A39). If & = 4, there is one way to connect each w to 
av, but this does not lead to a valid graph. 
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A37 A38 A39 


Case A.3.4 Assume u is not a 2-leaf. Then x1 is the only 2-leaf of H, so there is no x3. 
Then essentially the same argument as in Case A.3.2 repeats, with u and x, switching roles, 


and the same graphs are found. 


Case B. The stem set is S = {u,v,w}, and there are 2-leaves with neighborhoods {u, v}, 
{u,w}, and {v,w}. Thus u, v, and w will be colored black. 
Each 2-leaf of the inside graph H is in S, so H has at most three 2-leaves. 


Case B.1 If H has three 2-leaves, it may be K3 (B1). If not, it has order at least 4, so 
none of the 2-leaves of H are neighbors. Then each 2-leaf of G has distance more than 2 from 
a 2-leaf of H, which is impossible. 


Case B.2 If H has two 2-leaves, the third vertex in S' must be in both of their neighbor- 
hoods. Thus H has order at most 5. Thus H is kK, —e or Py + K, (B2, B3). 


Case B.3 If H has one 2-leaf v, then u must be one of its neighbors. If wu is a 2-leaf of 
H — v, H has order 5, so it is W;. There are two distinct choices for which vertex is w (B4, 
B5). If wu is not a 2-leaf of H — v, v has another neighbor, x, that is. Then u is adjacent to 
every vertex of H —v—v. If u is adjacent to x, then by Lemma 2.2, H is a 2-tree, so it has at 
least two 2-leaves, a contradiction. If u is not adjacent to x, then by Lemma 2.2, H — v — = is 
a 2-tree. 

Now z is adjacent to all 2-leaves of H —v—2, so H —v—z is a fan. Now w must be one 
of the 2-leaves of H —v— «x , but it cannot neighbor all vertices of the fan unless the fan is K3 
and H = W; , a previous case. 


Was EOS 


This completes the proof. 


A structural characterization of maximal 2-degenerate graphs with diameter 2 allows us to 
evaluate or bound parameters on this class, which would otherwise be difficult. Sharp bounds 
have been proved for the maximum degree of maximal planar graphs with diameter 2 [18, 20]. 
We state sharp bounds on the maximum degree A of maximal 2-degenerate graph with diameter 
2. A maximal 2-degenerate graph with A = n— 1 must have diameter at most 2. A maximal 
2-degenerate graph with A = n — 2 need not have diameter at most 2 (for example, add one 
vertex to a fan with at least 5 vertices). Proposition 1.2 implies 2-trees with diameter 2 have 
A> 2n, and this bound is sharp. 
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Corollary 2.6 A maximal 2-degenerate graph G with order n and diameter at most 2 has 


and this bound is sharp for all n. 


Proof For 1 <n < 4, there is only one maximal 2-degenerate graph, which has a dominating 
vertex. For n = 5, there are three such graphs, one (W; ) of which has no dominating vertex. 
The fact that maximal 2-degenerate graphs have size m = 2n—3 and minimum degree 2 implies 
A> 4 for n> 6. For 6 <n < 8, this is attained by adding 2-leaves to A4 and A23. 

Let G be a graph found under Case A, and H its inside graph. Then H has a stem that 
is adjacent to all 2-leaves of G with at most 5 vertices not adjacent to it, and only A39 attains 
this. Adding the 2-leaves of G to A39 as evenly as possible produces vertices with degree n — 6 
and n—-4— | 53]. Thus A > n-—6 for A39 when n > 12. Otherwise, A > n — 5 for graphs in 
Case A, and this is attained by graphs constructed from A37 when n > 9. 

Let G be a graph found under Case B, and H its inside graph with stem set {u,v, w}. 
Consider summing the degrees of u, v, and w. There are n — 3 other vertices, each of which 
is adjacent to at least two of u, v, and w. The graph induced by u, v, and w has at least two 
edges. Thus 2n — 2 = 2(n—3) +4 < d(u)+d(v)+d(w) < 3A, so A > [2 (n—1)]. This is 
attained by graphs constructed from B3. For n > 16, [2 (n—1)] <n —6, so the bound is as 
stated. 


We have seen that some maximal 2-degenerate graphs with diameter 3 are contained in a 


maximal 2-degenerate graph with diameter 2 (graphs A23-A39 above). The smallest maximal 
2-degenerate graphs not contained in a maximal 2-degenerate graph with diameter 2 have order 
7. They are all those with order 7 and diameter 3, excluding those listed in Theorem 2.5 (A25, 
A26, A29-A31, A33, A37, A38). 


Proposition 2.7 Let G be a maximal 2-degenerate graph. Then G is contained in a maximal 


2-degenerate graph with diameter at most 3. 


Proof If G has diameter at most 3, we are done. If not, consider a vertex v with maximum 
eccentricity. Let S be the set of all vertices with distance more than 2 from v. Add 2-leaves 
adjacent to v and each vertex in S, and call the set vertices added S’. Now the distance between 
v and any other vertex is at most 2. Vertices in S’ are all distance 2 from each other. A vertex 
in S’ and a vertex in G have distance at most 3, since there is now a path through v. Thus 


no new pairs with distance more than 3 are created. This process can be repeated with other 


vertices until a graph is constructed that contains G and has diameter at most 3. 
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§3. Diameter 2 k-Trees 


In this section, we prove a forbidden subgraph characterization of k-trees with diameter 2. 


Definition 3.1 A k-path graph G is an alternating sequence of distinct k- and k + 1-cliques 
€0, t1, €1,t2,°+* , tp, ep, starting and ending with a k-clique and such that t; contains exactly two 


k-cliques e;_1 and e;. 


Note that k-paths are also known a linear k-trees [1]. They are closely related to pathwidth 
[17]; in particular, they are the maximal graphs with proper pathwidth k. I have have further 
examined k-paths in two forthcoming papers [4, 5]. There is a simple characterization of k- 
paths. 


Theorem 3.2([15]) Let G be a k-tree withn > k+1 vertices. Then G is a k-path graph if and 
only if G has exactly two k-leaves. 


A k-path with a dominating vertex has nice structure. 


Lemma 3.3 A k-path has diameter at most 2 if and only if it has a dominating verter. When 
k > 2, ak-path with a dominating vertex can be represented as P+ K,, where P is ak—1-path. 


Proof Every k-path with order n < k +2 has diameter at most 2 and a dominating 
vertex. Consider constructing the k-path from K;, + K2, which has k-leaves u and v,, and 
N (u) = S; = N(v,). Iteratively add vertex v; with neighborhood $;, so that S; replaces one 
vertex of S;_; with v;_,. As long as $; and S; contain a common vertex, the graph has diameter 
2 and a dominating vertex. Once S; and $; do not contain a common vertex, the graph has 
diameter more than 2 and no dominating vertex. 

For the second claim, we use induction on order n. When n = k, G = Kx and the result 
holds. Let G be a k-path with order n > k containing a dominating vertex u, and assume the 
result holds for all graphs with order n—1. Then G has a k-leaf v, which is adjacent to u. Now 
G—v is a k-path with a dominating vertex, so it can be represented as P’ + K,, where P’ is 
a k —1-path. Then the other neighbors of v induce a clique in P’, so G can be represented as 
P+ kj. 


Note for k > 2, a k-tree with diameter 2 need not have a dominating vertex. 
Adding a k-leaf to a k-tree cannot change any existing distances. Thus when constructing 
a k-tree, the diameter can increase, but it cannot decrease, as it can in a maximal k-degenerate 


graph. 


Definition 3.4 A k-tree is minimal with respect to diameter 3 if deleting any k-leaf results in 


ak-tree with diameter 2. 


We can characterize these graphs. A tree is minimal with respect to diameter 3 if and only 
if it is Py. We have seen in Proposition 1.2 that a 2-tree is minimal with respect to diameter 3 
if and only if it is P?. In general, P3j,,5 is the smallest k-tree with diameter 3, but for k > 3 it 
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is not the only one. 


Algorithm 3.5 Let P be a k — 2-path, k > 3, of order n — 4 with k-leaves w and x. Join 
dominating vertices y and z to P, forming P+ K2. Add u with neighborhood Np (w) U {w, y}, 
and v with neighborhood Np (x) U {z, z}. Let Gz be the class of all graphs formed this way. 


z 


Theorem 3.6 A graph G is a k-tree minimal with respect to diameter 8 if and only if G € Gx. 


Proof (<=) Let G be a graph in Gx constructed using the algorithm. Then G is a k-tree, 
d(u,v) = 3, and u and v are the only pair with distance more than 2. 

(=) Let G be as stated. A k-tree with diameter 3 must contain a pair of vertices distance 
3 apart. Thus in a minimal k-tree with diameter 3, the vertices at distance 3 must be k-leaves, 
and no other vertices are k-leaves. Thus G is a k-path with leaves (say) u and v. Since G is 
minimal, G — u has diameter 2. By Lemma 3.3, it has a dominating vertex y, soG—u—yisa 
k — 1-path. Similarly, G— v has a dominating vertex z. Thus G — {u,v,y, z} is a k — 2-path. 
Then uw and v must each neighbor one of y and z, and one of the k-leaves of the & — 2-path. 


Thus G can be constructed using the algorithm, so G € Gx. 


Equivalently, a k-tree has diameter at most 2 if and only if it does not contain any graph 
in G,. When k = 3 and n > 8, the algorithm produces a unique 3-tree of order n minimal with 
respect to diameter 3 (shown below for n = 8). 
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Abstract: The total eccentricity of a graph G is defined as ¢(G) = > ecc(v) and the 
veV(G) 
Narumi-Katayama index of a graph G is defined as NK(G) = [][_ deg(v). In this paper, 
veV(G) 
we have obtained the total eccentricity and the bounds for Narumi-Katayama index for the 
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§1. Introduction 


Throughout this paper, all graphs we considered are simple and connected. For a vertex v € 
V(G), deg(v) denotes the degree of v. For vertices u,v € V(G), the distance d(u, v) is defined 
as the length of the shortest path between u and v in G. The eccentricity ¢(v) of a vertex v is 
the maximum among the distances from v to all the remaining vertices. The total eccentricity 
of the graph G, denoted by ¢(G) is defined as the sum of eccentricities of all the vertices of the 
graph G [3]. That is, 
¢(G) = S- ecc(v). 
vEV(G) 

In 1984, Narumi-Katayama [7] proposed a definition of a simple topological index which is 

defined as 
NK(G) = II deg(v). 
veV(G) 

On this graph invariant, several works [5,6,8-10] are reported and the name “Narumi-Katayama 
index” is used. 

In [8], I. Gutman et al. considered the problem of extremal Narumi-Katayama index 
and offered a few results filling the gap. For graphs without isolated vertices, I. Gutman et 
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al. [8] presented the minimal, second-minimal and third-minimal (maximal, second-maximal 
and third-maximal, resp.) NK-values for extremal graphs. Moreover, the maximal (second- 
maximal) Narumi-Katayama index of n-vertex tree (unicyclic graph) is determined [8] and the 
maximal Narumi-Katayama index of n-vertex bicyclic graphs is given. For connected n-vertex 
graphs, the minimal and second minimal Narumi-Katayama index are showed [8]. Consequently, 
the second-minimal Narumi-Katayama index among n-vertex trees and the minimal Narumi- 
Katayama index among n-vertex unicyclic graphs are presented [8]. 


KM. Kathiresan and 8. Arockiaraj introduced some generalization of complementary prism- 
s and studied the Wiener index of those generalized complementary prisms [9]. 


Let G and H be any two graphs on p, and pz vertices, respectively and let R and S be 
subsets of V(G) = {u1,u2,--- ,Up,} and V(H) = {v1, v2,--- , Up, } respectively. The comple- 
mentary product G(R)OH(S) has the vertex set {(ui,vj):1<i<pi,1< J < po} and (wi, v;) 
and (up, vz) are adjacent in G(R)OH(S) satisfying 


(i) ifi =h,u; € Rand vjv, € E(A), or ifi=h,u; ¢ R and vjy, ¢ E(#) or 
(it) if j =k,v; € S and uu, € E(G), or if j =k,v; ¢ S and ujun ¢ E(G). 


In other words, G(R)OH(S) is the graph formed by replacing each vertex u; € R of G 


by a copy of H, each vertex u; ¢ R of G by a copy of H, each vertex v; € S of H by a 
copy of G and each vertex v; ¢ S of H by a copy of G. If R = V(G) (respectively, S = 
V(#)), the complementary product can be written as GOH(S) (respectively, G(R)OH). The 
complementary prism GG obtained from G is GOK2(S) with |$| = 1. That is, GG has a copy 
of G and a copy of G with a matching between the corresponding vertices. 


In GG, we have an edge vd for each vertex v in G. The authors consider this edge as K or 
K1,1 or Py. By taking m copies of G and n copies of G, they generalize the complementary prism 
as a graph GOA(S), where H = Kmin (or Km n) and S is a subset of V(H) having m vertices 


and H = Com, (or Pam) whose vertex set is {v1,v2,°+* ,Vam} and S = {v1,0u3,+++ ,vam—1}- 


In [1,2], the eccentric connective index and first and second Zagreb indices have been bound 
for the generalized complementary prisms. 


Motivated by these works, the total eccentricity and the bounds for Narumi-Katayama 
index for the generalized complementary prisms of graphs have been found in this paper. 


§2. Main Results 


Proposition 2.1 Let G be any graph with k number of full degree vertices and k’ number of 
isolated vertices. Then for any m,n > 1, 
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Proof For any vertex v € V(Gm+4n), 


w) 2, if eccg(v) = 1 
ecc(v) = 
3, otherwise 


Hence, 


C(Gmtn) = ‘a ecc(v) 


vEV(Gmtn) 
m n 

= S- S- ecc(v) + :2 S- ecc(v) 
i=1 veEith copy of G t=1 veith copy of G 


= 38mp — mk + 3np — nk’ 


= 3p(m +n) — mk — nk’. 
Proposition 2.2 For any m,n > 1,¢(Gm.n) = 3mn. 


Proof For any vertex v € V(Gm,n), its eccentricity is 3. Hence 


¢(Gm,n) = Se ecc(v) = 3mn. 


vEV(Gin,n) 


Proposition 2.3 For any m > 2,C(Gh, m) = 2mp(m + 1). 


Proof In Gf, m,ecc(v) =m + 1, for all v € V(G§, ,,). Hence, 


m,m? 


C(Gom)= Sd. ece(v) = 2mp(m +1). 


vEV(GS, m) 
Proposition 2.4 For anym > 1, ¢(GI, ,) = mp(3m + 1). 


Proof Let v;;,0j,7,1 <j < p be the vertices of the i” copy of G and G respectively in 
GF nm for 1 <i <m. Notice that for 1 <i< [#] and1 <j <p, ecc(uij) = 2m +2 — 2%; for 
1<i< |@| andl <j <p, ecc(vjj) = 2m+1- 2%; for [Z]+1<i< mand1<j <p, 
ecc(vj,;) = ecc(Um+41-i,j); for |B] +1<i< mand 1<j <p, ecc(tjj) = ecc(vm41-ij)- 


From these, 


(Grim) = >) ece(v) 
vEV(GE, om) 
= 2p[(m +1) + (m+ 2) +--+ + 2m] 
= Qp[m? +14+24+---+m] 
m(m + 1) 
2 


= 2p rn? + | = mp[3m + 1]. 


Proposition 2.5 If G,G € Fr, then ¢(GG) = 4p. 
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Proof If G and G are members of F22, then any two vertices with in G or with in G are 
with in distance 2. If u and v are at distance 2 in G, then U 0, ut € E(GG) and hence u and U 
are at a distance 2 in GG. Hence ecc(v) = 2 for all v € V(GG). Therefore, 


¢(GG) = a ecc(v) = 4p. 


veEV(GG) 
Proposition 2.6 If G € F,,, then ¢(GG) = 5p. 


Proof If G € F,;, then G is a complete graph and G is a totally disconnected graph. So 
GG is simply K, 0 Ky. 


In GG, the vertices on the copy of G are of eccentricity 2 and vertices on the copy of G 


are of eccentricity 3. Hence ¢(GG) = 5p. 


Proposition 2.7 If G € Fiz, then ¢(GG) = 5p. 


Proof Let u be the full degree vertex in G. If v is any vertex other than u in G, then 
d(u,v) = 1 and d(u,d) = 2. So ecc(u) = 2 in GG. As 7 is an isolated vertex in G, by the edges 
wu, wv and vv, any vertex U in the copy of G will be reached at a distance 3 from @ and any 
vertex v in the copy of G will be reached at a distance 2. So ecc(v) = 3 in GG for all 0 in the 
copy of G. 

Let v be a vertex of eccentricity 2 in G. Whenever z is a vertex of distance 2 from v in G, 
the edges «% and Z U give that d(v,Z) = 2. Whenever x is an adjacent vertex to v, the edges 
vax and x% give again that d(v,Z) = 2. Hence ecc(v) = 2 in GG. So for any vertex v € V(GG), 


2, if v is in the copy of G 


ecc(v) = i 
3, if v is in the copy of G. 


Hence, 


¢(GG) = S- ecc(v) = 2p + 3p = 5p. 
veEV(GG) 


Proposition 2.8 If G € F53, then ¢(GG) = 5p. 


Proof If v is a vertex of eccentricity 2 (or 3) in G, then v is also of eccentricity 2 (or 3) in 
GG. But in G, the vertex 0 corresponding to v in G is of eccentricity 3 (or 2) and hence 7 is 
also of eccentricity 3 (or 2) in GG. That is, 


2: if eceg(v) = 2 or eccg(v) = 
ecc(v) = 


3 
3, if eccg(v) = 3 or ecca(v) =3 
So exactly p number of vertices are of eccentricity 2 in GG and the remaining p vertices are of 
eccentricity 3 in GG. Hence ¢(GG) = 5p. 


Proposition 2.9 If G € Fo4 with k number of vertices of eccentricity 2, then ¢(GG) = 6p — k. 
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Proof When G € Fo4, G € Fo2 and hence ecc(v) = 2 in GG whenever u and v are at 
distance more than 2 in G, by the edges u U,U 0, dv in GG, d(u,v) = 3 in GG. Also if u is a 
vertex of eccentricity 2 in G, then it is also u is a vertex of eccentricity 2 in GG. Thus exactly 
k number of vertices are of eccentricity 2 in GG and the remaining vertices are of eccentricity 
3 in GG. Hence ¢(GG) = 2k + 3(2p — k) = 6p — k. 


Proposition 2.10 If G € F3, then ¢(GG) = 5p. 


Proof If G € F3, then G € F2 and hence ecc() = 2 in GG whenever T is in the copy of 
G. Also since eccg(v) > 3, ecc(v) = 3 in GG. Hence ¢(GG) = 5p. 


Proposition 2.11 For any positive integers m and n, 


degg(u) +m+n-—1, when w is in a copy of G 


dega(u) +m+n-—1, when wu is in a copy of G. 


We know that A.M. > G.M. and a+b > 2Vab. Therefore 


NKGmin)=  T] deg(u) 
uEV(Gm+n) 
= II deg(u). II deg(t) 
u€ copies of G WE copies of G 
= II deg(u) | . Il II deg(u) 
i=1 \ue ith copy of G t=1 \ae ith copy of G 
= Il (deg(u) +m+n-—1) II (m+n—1+ dega(t)) 
ueV(G) ueV(G) 
> II 2\/dega(u)(m +n — 1) II 2y/(m-+n- 1)degg(t) 
ueV(G) ueV(G) 


m 


2 


=(2Vm4+n—1)""{ J] dega(u) | (2¥m+n—1)"{ [J dege(a) 


ueV(G) ucV(G) 


Ns 


= [2Vm+n—1\°"*”” INK(G)|? [NK(G@)]3. 


Proposition 2.12 For any positive integers m and n, 


n 
m nm 
2 2 


NK(Ginn) > 2'°"*™? (nF m?)? (NK(G))? (NK(G))?. 
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Proof In the graph Gmn, 


di 
deg(u) = ega(u) 
dega(u) 
Therefore, 
NK(Gmn)= [I deg(u) 


uEV(Gm,n) 
ue copies of G 
m 


i=lue it? copy of G 


deg(u). II 


+n, when u is in a copy of G 


+m, when wu is in a copy of G. 


deg(t) 


UWE copies of G 


deg(u)| . II II 


t=1Ge ith copy of G 


deg(t) 


=| TT Gegatu)+n)} .| TI] ea) +m) 
ueV(G) ueV(G) 

>| TI evdesatan) | I] @yaese(a)m) 
ueV(G) meV (G) 


(2vn)"" | [J dege(u) 


ueV(G) uEeV(G) 
= 2("+")P (nF m2)” NK(G)?NK(G)?. 
Proposition 2.13 For any positive integers m > 2, 
NK(Giim) = 8?[NK(G)NK(G)]?. 


Proof In the graph GS 


mim? 


Ns 


deg(u) 


dega(u) + 2, when uw is in a copy of G 
deg(u) = ee 
dega(u) + 2, when uw is in a copy of G 
Therefore, 
NKGrm)= TT adgm= JT deg. I] 
ueEV (GS, an) ue copies of G mE copies of G 


-T Il 


ue it? copy of G 


deg(u)| . II II 


ueith copy of G 


deg(u) 
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- II (dega(u) +2)] . II (dega(w) + 2) 
ueV(G) ueV(G) 

>} [] @v2dege(u))| -| [] @y/2degg@)) 
ueV(G) meV (G) 


II 
ae 
SI 
s 
= 
Q 
g 
Q 
Ss 
time 
Ss 
Se 
3 
8 
a 
Q 
g 
QI 
a 


This completes the proof. 


Proposition 2.14 For any positive integers m > 2, 


mm 
2 


NK(Giim) = 2? (NK(G@)NK(G)] 


Proof In the graph G? 


mim? 


degg(u) +1, when uw is in the first copy of G 
degg(u) + 2, when uw is in the remaining copies of G 
deg(u) = oe 
dega(u) +1, when uw is in the last copy of G 
dega(u) + 2, when wu is in the remaining copies of G 
Therefore, 
NK(Ginm)= [[ — deg(u) 
ueV(GE, om) 
= II deg(u). Il deg(%) 
ue copies of G UE copies of G 
= II II deg(u) } . II II deg(t) 
i=lue it? copy of G i=1 Ge ith copy of G 
~ TL deo | (TE = TL aeatwy 
ue first copy of G i=1ue ith copy ofG 


T TL ccs TL cata) 


t=1 Geith copy ofG ue mt copy of G 
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= Il (dega(u)+1)] | [J [] (ega(u) + 2) 


ue first copy ofG i=2 ueV(G) 
i. (dega(u@) + 2) |] . II (dega(u) + 1) 
i=1 G) ue mth copy ofG 
m-1 
>| J] 2Vdega(u) I] 2v2dega(u) 
ueEV(G) ueEV(G) 
m-1 
Il eee Tw) Il 2,/dega(t 
uEeV(G ueV(G) 
= ae Tr (dega(u))? [[ (dega(u)) "= 
ueV(G) ueEV(G) 
II (dega(u aes (dega(u ))2 
ueV(G) EV(G) 
- 7 


Sgr? II deg(u) II deg(t) 


ueV(G) uceV(G) 


= 2?""! |NK(G)NK(G)| 


mm 
2 
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